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ABSTRACT 

Recent studies have shown that various inaccuracies occur when 
a single -degree -of-freedom gyroscope is subjected to an angular vibra- 
tory environment. Under such conditions the performance of the gyro- 
scope has been analyzed and the nature of the inaccuracies investigated. 
Erro'rs due to "coning 11 , ’’dynamic unbalance 11 , and "sculling 11 have been 
included in the investigation. 

Values of the inaccuracies representative of those expected from 
tests using the precision angular vibrator developedby the Instrumentation 
Laboratory of the Massachusetts Institute of Technology are presented. 

By reference to a dynamic model the equation of motipn was ob- 
tained. A suitable solution to the equation was determined in series 
form for the case of sinusoidal inputs simultaneously applied about two 
axes. No restrictions are placed on the frequency of the applied mo- 
tions. The terms in the series solution reveal the nature of the angular 
displacement of the float with respect to inertial space about the output 
axis for two cases : (1) zero elastic restraint, (2) finite elastic restraint. 

It is strongly recommended that correlation of actual experimen- 
tal results with the developed theory be obt ained prior to further 
theoretical considerations. 



Thesis Supervisor: Dr. Walter Wrigley 

Title: Professor of Instrumentation and 
Astronautics 



li 



ACKNOWLEDGEMENT 



The authors wish to express their appreciation to the 
personnel of the Instrumentation Laboratory, Massachusetts 
Institute of Technology, who assisted in the preparation of 
this thesis. Particular thanks are due to Dr. Walter Wrigley, 
the thesis supervisor; Mr. Arnold J. Smith, who as technical 
advisor, guided and inspired the entire project ; Mr. Herbert 
A. Weinstock, who furnished considerable technical assistance; 
Mr. Daniel Goldenberg, for his mathematical contributions; and 
Mr. Peter J. Palmer, who suggested the thesis topic, for his 
cooperation. 

Acknowledgement is more than due to Miss Nancy Keene 
who typed and corrected the original manuscript. 

The graduate work for which this thesis is a partial require- 
ment was performed while the authors were assigned to the U. S. 
Naval Administrative Unit, Massachusetts Institute of Technology. 

This report was prepared under Project 53-173, Division of 
Sponsored Research M. I. T. , sponsored by the U. S. Navy Special 
Projects Office under Bureau of Ships Contract NObs 78136. 



iii 



TABLE OF CONTENTS 



Chapter Page 



Object 1 

1 INTRODUCTION 2 

2 THE FLOATED SINGLE -DEGREE -OF- 

FREEDOM INTEGRATING GYROSCOPE. ... 4 

2. 1 Design Features 4 

2. 2 Operation 5 

3 INACCURACIES DUE TO SIMULTANEOUS 
VIBRATIONS ABOUT THE CASE INPUT 

AND SPIN AXES 12 

3. 1 General 12 

3. 2 Coordinate System and Symbols 12 

3. 2. 1 Coordinate System 12 

3. 2. 2 Definition of Symbols 14 

3.3 Gyro Response to Angular Vibrations .... 16 

3. 3. 1 Applied Torque Due to Angular 

Vibrations 16 

3. 3. 2 Dynamic Torque 17 

3. 3. 3 Gyroscopic Torque 19 

3. 4 Coning 19 

3. 4. 1 General Discussion of Coning 19 

3. 4. 2 Platform Coning Effect 20 



IV 



3. 5 Sculling 26 

4 MAGNITUDE OF THE ANGULAR VIBRATION 

INACCURACIES 29 

4. 1 General 29 

4.2 Magnitude Determination 30 

4. 2. 1 Magnitude of the Dynamic Torque . . 30 

4. 2. 2 Magnitude of the Gyroscopic Torque . 32 

4. 2. 3 Magnitude of Platform Coning ... 36 

4. 2. 4 Magnitude of Sculling Effect ... . 37 

5 THE RESPONSE OF THE SINGLE -DEGREE - 

OF-FREEDOM INTEGRATING GYROSCOPE 
TO TWO ANGULAR VIBRATIONS SIMUL- 
TANEOUSLY APPLIED ABOUT THE SPIN 
REFERENCE AND OUTPUT AXES .... 41 

5. 1 Introduction 41 

5. 2 Definition of Symbols 42 

5. 3 Equation of Motion of the Dynamic Model . . 43 

5. 4 Method of Solution of the Equation of Motion 
for Two Different Sinusoidal Vibration 
Inputs of Small Magnitude 44 

5. 4. 1 Case I; Negligible Elastic Restraint 

Between the Float and Case About 

the Output Axis 46 

5. 4. 2 Case II; Finite Elastic Restraint 
Between the Float and Case About 
the Output Axis 50 

5. 5 Drift of a Typical Gyroscope for the Case of 
Negligible Elastic Restraint About the 
Output Axis with Two Equal Input 

Frequencies 53 



v 



6 CONCLUSIONS AND RECOMMENDATIONS. . 57 

6. 1 Conclusion 57 

6. 2 Recommendations 60 

Appendix 

A SERIES SOLUTION FOR THE EQUATION OF 

MOTION OF A SINGLE -DEGREE -OF- 
FREEDOM INTEGRATING GYROSCOPE . . 62 

A. 1 Series Solution for Negligible Elastic 
Restraint Between the Float and Case 
About the Output Axis ( k 0 ) 62 

A. 2 Series Solution for Finite Elastic Restraint 
Between the Float and Case About the 
Output Axis (k 0) 75 



vi 



LIST OF ILLUSTRATIONS 



Figure Page 

2-1 Pictorial Schematic of the Single -Degree -of 

Freedom Floated Integrating Gyro Unit ... 6 

2-2 Sectioned View of a Typical Floated Single - 

Degree-of-Freedom Integrating Gyro ... 7 

2-3 Line Schematic Diagram for the Single - 

Degree-of-Freedom Integrating Gyro Unit . . 8 

2- 4 Functional Diagram for the Floated 

Single -Degree -of- Freedom Integrating 

Gyro Unit 9 

3- 1 Gyroscope Coordinate System and 

Angular Velocity Components 13 

3- 2 Euler Angle Relations Between Platform 

and Desired Orientation Coordinates .... 21 

4- 1 Drift Rate Due to Cross Coupling Coning for the 

Single-Degree-of-Freedom Integrating Gyro . 34 

4-2 Estimated Platform Coning Drift 38 

4- 3 Estimated Sculling Drift 40 

5- 1 A Dynamic Model of the SDFI Gyroscope for 

Motions About the Spin and Output Axes ... 45 

5-2 Drift Rate vs. Frequency for a Typical Gyro . 56 



OBJECT 



The object of this thesis is to analyze the performance 
of a floated single -degree -of-freedom integrating gyroscope 
subjected to angular vibrations in order to determine the 
inaccuracies that may result from such an environment. 



CHAPTER 1 



INTRODUCTION 



The age of the golden sixties is upon us. Ushered in with this 
decade is an era of accomplishment which has raised the importance 
of the gyroscope to an all-time high. This era of great advancement 
includes those fields of Science and Engineering which are brought to 
the fore when one thinks of space. Because of space, these fields 
have been searched diligently for methods of navigation and control 
and from this searching the gyroscope has emerged, holding a posi- 
tion of extreme prominence. 

With this relatively new stature, the gyroscope has naturally 
been investigated in efforts to improve its performance as well as 
to predict what may be expected in the future. The floated single - 
degree-of-freedom integrating gyroscope, an important element in 
most space navigation and guidance systems, is the subject of this 
analysis. The SDFI gyroscope has been interrogated from the 
standpoint of being subjected to varying vibratory inputs which may 
introduce error. These vibrations may be of very low frequency 
due to the bending modes of long missile bodies or they may be of 
high frequencies associated with aerodynamic flutters and buzzing. 
Through the theoretical results, it is hoped that undesirable effects 
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caused by such vibratory motions will be recognized and later minimized 
when actual verifying data is available. 

The floated single -degree -of -freedom integrating gyroscope is 
an instrument capable of sensing angular velocities to a high degree 
of accuracy. When utilized in inertial guidance systems, the gyro- 
scope must have considerable stability and extreme sensitivity. The 
vibratory environment in which the gyro functions may introduce un- 
wanted angular outputs. This degradation of desired performance is 
the result of both linear and angular oscillations. An analysis of the 
undesired effects is made in order to determine the nature and magni- 
tudes of the errors involved. 

This investigation was conducted with the view that the Massachu- 
setts Institute of Technology Instrumentation Laboratory will conduct 
tests on a precision angular vibrator. A previous investigation has 
disclosed what may be expected if two equal frequencies are simul- 
taneously applied to the output and spin reference axes of a typical 
single -degree -of -freedom gyroscope. The present analysis is concerned 
with any two frequencies simultaneously applied to the output and spin 
reference axes of a SDFI gyroscope. 

The existing equation of motion is based on a simplified dynamic 
model which accounts for damping, compliance, and inertia. As no 
closed form solution to the equation is known, a series form was used 
which gives a very clear representation of the gyroscopic drift terms. 

The series solution appears easy to handle and analyze in view of the 
theory of superposition. The terms seem reasonable in magnitude and 
should be easily verified by the Massachusetts Institute of Technology 
angular vibrator. 



3 



CHAPTER 2 



THE FLOATED SINGLE -DEGREE -OF-FREEDOM 
INTEGRATING GYROSCOPE 



2. 1 Design Features 

The types and locations of components within the ideal floated 
single -degree -of -freedom integrating gyroscope unit is pictorially 
shown in Fig. 2-1. These components are the gyroscopic element, 
gimbal, viscous damper, microsyn torque generator, microsyn signal 
generator, gimbal bearings, and the case. 

The gyroscopic element and gimbal are hermetically sealed in 
a cylindrical shell which is separated from the outer case by a high 
density fluid. The case is also hermetically sealed. The two shells 
form a thin-walled cylinder that performs two vital functions when 
completely filled with a high density fluid: 

1. it provides a viscous damping torque that opposes the 
precessional torque generated by the gyro element, 

2. it provides a buoyant force that neutralizes the weight 
of the gyro element and its supporting structure. 

The clearance between the gimbal float and the case plus the fluid 
viscosity determines the damping of the viscous damper; the fluid 
density and the amount of displaced fluid determine the buoyant force. 
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Power is conducted to the hysteresis type alternating -cur rent 
motor by flex leads. The entire float assembly, or torque summing 
member, is initially aligned with the case by jewel bearings. Final 
alignment is by magnetic suspension. The suspension (elastic re- 
straint to float radial motions) is provided by the microsyns in con- 
junction with their primary functions of angle measurement and torque 
generation. The magnetic forces necessary for suspension originate 
with the change in excitation currents of properly tuned microsyn 
coils as a result of the radial displacement of the microsyn rotor from 
its centered position. Figure 2-2 is a sectioned view of a typical 
floated single-degree-of-freedom integrating gyroscope. A more de- 
tailed discussion of the single-degree-of-freedom integrating gyroscope 
may be found in Refs. 1 and 2. 

2. 2 Operation 

An ideal single-degree-of-freedom integrating gyroscope receives 
as inputs, the angular velocity of its case with respect to inertial space 
about the input axis and electrical current input signals to the torque gen- 
erator. The unit generates an output voltage proportional to the inputs. 

In open loop operation, the gyro output signal would represent the angu- 
lar deviation of the gyro case from a reference orientation. The refer- 
ence is established from initial conditions and by command signals to 
the torque generator. 

In operation, the gyro unit is housed in a thermally controlled 
environment that maintains the temperature to within a few degrees 
Fahrenheit. This control is important because temperature affects the 
density and viscosity of the fluid. 
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Fig. 2-1. Pictorial schematic of the single-degree-of-freedom floated integrating gyro unit. 



Float-Girobal Structure Heater and Thermostat 

Gyro Mount \ Sensing Element 




7 



Fig. 2-2 Sectioned view of a typical floated aingle-degree-ef-freedem 

integrating gyro. 




CASE, (co) or (c) 

THE STRUCTURE THAT GIVES-SUPPORT FOR THE 
INTERNAL WORKING PARTS OF THE GYRO UNIT, 
ENCLOSES THESE PARTS, AND CARRIES PROVI- 
SIONS FOR EXTERNAL CONNECTIONS OF ALL KINDS. 

TORQUE GENERATOR AND MAGNETIC SUSPENSION 
UNIT, (tg) + (ms) 

COMPONENT FOR (a). RECEIVING INPUT SIGNALS AND 
PRODUCING A CORRESPONDING OUTPUT TORQUE AP- 
PLIED TO THE GIMBAL ABOUT THE OUTPUT AXIS AND 
(b) ALSO PROVIDING MAGNETIC FORCE ACTING TO 
CENTRALIZE THE PIVOT SHAFT IN THE BEARING. 
DAMPER, (dmp) or (d) 

SUBSYSTEM RECEIVING THE ANGULAR VELOCITY OF 
THE GIMBAL WITH RESPECT TO THE CASE AS ITS INPUT 
AND PRODUCING AS THE OUTPUT A RETARDING TORQUE 
ACTING ON THE GIMBAL ABOUT THE OUTPUT AXIS WITH 
A MAGNITUDE PROPORTIONAL TO THE MAGNITUDE OF 
ITS ANGULAR VELOCITY INPUT. 



GYRO UNIT, (gu) 

THE ENTITY MADE UP OF THE COMPONENTS REPRESENTED 
IN THIS DIAGRAM AND ALL THE ADDITIONAL PARTS NECES- 
SARY FOR A SINGLE PACKAGE TO CARRY OUT THE FUNC- 
TIONS OF A GYRO UNIT. 

SIGNAL GENERATOR AND MAGNETIC SUSPENSION UNIT, 
(sg) + (ms) 

COMPONENT FOR (a) RECEIVING THE ANGLE OF THE SPIN 
AXIS WITH RESPECT TO THE CASE AS ITS INPUT AND PRO- 
DUCING A CORRESPONDING SIGNAL THAT SERVES AS THE 
OUTPUT SIGNAL FROM THE GYRO UNIT AND (b) ALSO PRO- 
VIDING MAGNETIC FORCE ACTING TO CENTRALIZE THE 
PIVOT SHAFT IN THE BEARING. 

GIMBAL, (gim) 

STRUCTURE CARRYING THE BEARINGS FOR THE SPINNING 
ROTOR OF THE GYROSCOPIC ELEMENT, ROTORS FOR THE; 
TORQUE GENERATOR AND SIGNAL GENERATOR, PART OF 
THE DAMPER, FLOAT SEALS AND STRUCTURE, BALANCE 
ADJUSTMENTS, STOPS, PIVOTS, ETC. 

RIGID ARM UNBALANCE COMP ENSATOR.(rauc) 

SCREW ADJUSTMENTS FOR SHIFTING THE CENTER OF MASS 
OF THE GIMBAL FLOAT TO A POSITION SUBSTANTIALLY ON 
THE OUTPUT AXIS. 



* A DISCUSSION OF GENERALIZED CONVENTIONS FOR SELF-DEFINING SYMBOLS OF WHICH Aj (r<jf IS AN EXAMPLE IS GIVEN BY 
DRAPER, McKAY AND LEES IN INSTRUMENT ENGINEERING (10), VOL. I. 

Fig. 2-3 Line schematic diagram for the single-degree-of-freed*m 
integrating gyra unit. 
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Fig. 2-4 Functional diagram for the floated single-degree-of-freedoa 
integrating gyro unit. 



The torque and signal generators are microsyn units. The rotors 
of these units are attached to the gimbal float while the stators are fixed 
to the case. 

A line schematic diagram is shown in Fig. 2-3 and a functional 
diagram in Fig. 2-4. If the gyroscopic element output is the only torque 
applied to the torque -summing member of the gyro unit, it is balanced 
by the viscous shear torque. The viscous torque is proportional to the 
angular velocity of the floated gimbal assembly relative to the case. 

This torque balance causes the output angular velocity to be proportional 
to the input angular velocity. Thus, the angle of the floated gimbal with 
respect to the case is proportional to the angular displacement of the 
case about the input axis, measured with respect to inertial space. In 
other words, the angle through which the float turns with respect to the 
case represents the time integral of the gyro input angular velocity. 

The torque generator permits introduction of command orienta- 
tion signals. The torque generator applies a torque to the float, propor- 
tional to the current supplied the torque generator. If there is no 
gyroscopic torque, the output of the torque generator causes the float 
to rotate (against the viscous shear torque of the damper) at a rate 
proportional to the current input to the torque generator. The gyro 
output is therefore equivalent to an angular displacement about the gyro 
input axis. The output axis angular displacement is measured by the 
microsyn signal generator. The output voltage of the signal generator 
is utilized for indication or control. 

For the gyroscope to be a linear instrument, the angle through 
which the float turns with respect to the case must be kept small. For 
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the integrating gyro this implies closed loop operation. There are 
two types of closed (servo) loops presently in common use: 

1. a rate feedback loop in which the signal generator output 
voltage is amplified and applied to the gyro torque gener- 
ator to act as an elastic restraint, 

2. a stabilization servo loop in which the signal generator 
output voltage is amplified and applied to a servo motor 
that supplies an angular rate about the gyro input axis. 
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CHAPTER 3 



INACCURACIES DUE TO SIMULTANEOUS 
VIBRATIONS ABOUT THE CASE INPUT AND SPIN AXES 

3. 1 General 

This chapter deals with the response of the single-degree-of- 
freedom integrating gyroscope when simultaneous angular vibrations 
are applied about the case input and spin axes. The resulting error 
torques are determined and their physical origin analyzed. 

It is assumed that when an angular vibration is present in the 
plane of the gyro case input and spin axes, components exist about both 
axes. The torque applied externally to the gyroscope is developed using 
basic gyro theory. 

3. 2 Coordinate System and Symbols 
3. 2. 1 Coordinate System 

The coordinate systems and the relationship of the angular ve- 
locity components are shown in Fig. 3-1. The spin reference axis 
(SRA), input reference axis (IRA) and output reference axis (ORA) are 
orthogonal axes fixed to the gyro case. The gyroscopic element is 
considered to be rigidly fixed to the float so that the float input axis (IA) , 
spin axis (SA) and output axis (OA) are also orthogonal. 
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OA 




W (IF)SA " W (IC)SA cos A (gim) + W (IC)IA Sln A (gim) 
W (IF)IA = W (IC)IA cos A (gim) “ W (IC)SA sln A (gim) 



Fig. 3-1. Gyroscope Coordinate System and Angular Velocity 

Components 
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3. 2. 2 Definition of Symbols 



The following symbols and definitions will be used in the analysis. 



A (C-F)OA" A (gim) 


Angle of float with respect to case 


IA 


Float input axis 


IRA 


Case input axis 


ORA = OA 


Output axis of both the case and float 


SA 


Spin axis 


SRA 


Spin reference axis (coincident with SA 
when A (gim) = 0) 


W SP 


Angular velocity of gyro wheel about 
spin axis 


W (IC)IA 


Angular velocity of case with respect to 
inertial space about case input axis 


W (IF)IA 


Angular velocity of float with respect to 
inertial space about float input axis 


W (IC)SA 


Angular velocity of case with respect to 
inertial space about case spin axis 


W (IF)SA 


Angular velocity of float with respect to 
inertial space about float spin axis 


J IA 


Moment of inertia of wheel and float 
about float input axis 


J SA 


Moment of inertia of wheel and float 
about float spin axis 
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H 

H NS 

H SP 



M (app) 



dH 



dt 



dH 

dt 



F 



M (out) 

i. J> k 



n 



n 



x 



n 



z 






Moment of inertia of gyro wheel about 
float spin axis 

Total angular momentum 

Nonspin angular momentum 

Angular momentum of spinning 
wheel 

Torque applied to the gyro unit 
input axis 

Rate of change of angular momentum 
with respect to inertial space 

Rate of change of angular momentum 
in float coordinates 

Torque output about gyro output axis 

Unit vectors along orthogonal axes 
forming right-hand set 

Linear acceleration 

Linear acceleration component along 
platform axis having i as unit vector 

Linear acceleration component along 
platform axis having j as unit vector 

Linear acceleration component along 
platform axis having k as unit vector 

Linear acceleration components along 
non-rotating axes having i, j, k as unit 
vectors 

Euler angles 



15 



3.3 



Gyro Response to Angular Vibrations 



3. 3. 1 Applied Torque due to Angular Vibrations 

In vector notation the torque applied^ to the gyro is 



(app) 



dH 




dH 


dt 


I 


dt 



+ W (IF) X H 



(3-1) 



By Newton's Law of Action and Reaction, the gyro unit output torque is 



M (out) = - M (app) = H x W (IF) 




(3-2) 



Using the relationships shown in Fig. 3-1, the following vectors may be 
written: 

W (IF) = 1 W (IF)SA + j W (IF)IA 

H = H NS + H SP = 1 ^SA W (IF)SA + X W W SP^ + ^^A W (IF)D?) 

W (IF) X H = k ^ I IA _I SA ) W (IF)SA W (IF)IA + H SP W (IF)LA 1 (3-3) 

The relationships from Fig. 3-1 are then substituted into 
Eq. (3-3) to obtain Eq. (3-4). 

M (app) = k t^IA ~ I SA ) ( cos2 A (gim) " sin2 A (gim) ) (W (IC)SA W (IC)IA? 

+ ^LA " W ^ w 2 (IC)IA " w2 (IC)SA^ sln A (gim) cos A (gim) 

- Hg p (W( IC)IA cosA^ gim) - W^ IC)SA sinA^ g . m) ) ] (3-4) 

For a practical single-degree-of -freedom gyroscope the gimbal angle, 
A(gim) usually does not exceed one degree. v ’ Therefore, the equation 
may be simplified by making use of small angle approximations and 
neglecting second order terms in The torque equation is 

reduced to: 
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X 



M (app) " k I (I I A I SA ) W (IC)S A W (IC)IA + (I I A - W A (gim) 

(W (IC)IA " W (IC)SA ) ” H SP (W (IC)IA ~ A (gim) W (IC)SA> 1 

(3-5) 

For a single -degree -of-freedom gyroscope only torques about 
the output axis are considered. The rate of change of angular momentum 
in float coordinates does not contribute an output axis torque. Equation 
(3-5) represents the applied torque of the gyro and, since it is all along 
a single axis, vector notation will not be necessary. 

M (app) = t (I IA " W W (IC)S A W (IC)I A + (I IA ~ W A (gim) (W (IC)IA ~ W (IC)S^ 
- Hg p (W^ IC)IA - A (gim) (W( IC)SA ) ] (3-5a) 

The terms of Eq. (3-5) represent error or undesired torques. 

A single -degree-of-freedom gyroscope, operating in a closed loop, 
would have an angular displacement from the null or zero position that 
would be in error due to these undesired torques. An examination of 
these error torques will now be made. 

3. 3. 2 Dynamic Torques 

The first torque component to be considered will be 

°IA-W W (lOSA W aC)IA (3 - 6) 

It is evident that the nature of the error will depend upon the variation 

of the two simultaneous angular velocities and the float inertias about 

f5i 

the spin and input axes. This error is known as anisoinertia v ' or 
dynamic unbalance. It could be eliminated by building gyro units 
having equal inertias about both axes. 
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Next consider the second term of Eq. (3-5) 

^IA “ W A (gim) ^ W (IC)IA “ W (IC)SA^ ^ 3 ‘ 7 ^ 

In the above expression the nature of the error will be determined by 
the difference of the inertias as well as the difference of squared angu- 
lar velocities. Variation of gimbal angle must also be considered since 
it is not a constant. 

Expressions (3-6) and (3-7) are dynamic torques^ and may be 
interpreted as representing four angular momentums and their associa- 
ted precessional angular velocities. Consider that there are four angular 
velocity components and four corresponding angular momentums. 



Momentum Components 


Angular Velocity Components 


X SA W (IC)SA 


W (IC)SA 


" A (gim) J IA W (IC)SA 


" A (gim) W (IC)SA 


I IA W (IC)IA 


W (IC)IA 


A (gim) X SA W (IC)IA 


A (gim) W (IC)IA 



The total ’’dynamical torque”^ when each forced precession is 
considered is: 

M D = “ I SA W (IC)SA (W (IC)IA~ A (gim) W (IC)SA^ 

" I IA A (gim) W (IC)SA ( W (IC)SA + A (gim) W (IC)IA ) 

+ J IA W (IC)IA (W (IC)SA + A (gim) W (IC)Ia) 

“ J SA A (gim) W (IC)IA (W (IC)IA “ A (gim) W (IC)SA^ 
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The part in parentheses for each of the above terms is the forced pre- 
cessional rate while the part before the parentheses is the angular 
momentum. For example, the angular momentum component, IgA^(IC)SA 
precesses about the input axis with the angular velocity, 

^(IC)IA ~ ^(gim)^(IC)SA° When the equation is expanded and second 



order terms of A, 



(gim) 



neglected, the expression for dynamic torque is 



obtained. 

M D = (I IA~ W W (IC)IA W (IC)SA + (I IA~ I SA ) A (gim/ W (IC)IA~ W (IC)SA^ 

(3-8) 

3. 3. 3 Gyroscopic Torque 



The applied gyroscopic torque on the gyro is 
M ge = ' H SP (W (IC)IA ' A (gim) W (IC)SA ) = ' H SP W (IF)IA 

(3-9) 

Variation of this torque will, as in previous considerations, be de- 
pendent on the angular velocities and the gimbal angle. The component, 
Hgp A (gi m ) w (ic)SA 9 is s P ecial significance and is called Coning 
Error. ^ 

3. 4 Coning 

3. 4. 1 General Discussion of Coning 

The coning error defined above is more usually considered a 
cross coupling error. It is a result of oscillations in the plane of the 
case input and spin axes. Comparing the form of the coning expression 
from the next section with that of cross coupling reveals a close similar- 
ity. Because of this similarity, it is possible to consider cross coupling 
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as another cause of the coning error. This form of coning error may 
be reduced by restraining rotations about the output axis. A tight servo 
loop would keep the error negligible. 

If the gyro input axis has a motion such that it returns periodically 
to its original position, then the input axis will trace out a path in space. 
This motion is known as the coning motion and is present even in a theo- 
retically perfect instrument. 

The coning error is a geometric effect that allows the gyro to 

(7) 

respond to an input rate produced by the coning motion. Almost any 
motion of the case creates coning. Consider oscillations in the plane 
of the case input and spin axes. The component along the input axis 
will cause the gyro to precess about the output axis. When the oscillating 
motion about the spin reference axis is combined with this precession, the 
result is a coning motion of the float input axis. As another example, if 
there is oscillation of the case in the spin axis - output axis plane, coning 
would again result because the damping fluid will drag the float around. 

When single-degree-of-freedom gyros are used for stabilization, 
a coning error exists that cannot be eliminated. 

3. 4. 2 Platform Coning Effect ® 

Consider two coordinate systems: 




platform coordinate system 




coordinate system defining the 
orientation the platform servos 
are trying to maintain 
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Fig. 3-2. Euler Angie Relations Between Platform and Desired 

Orientation Coordinates 
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By use of the Euler angle rotations^taken in the order shown in 
Fig, 3-2, the coordinate transformations, using matrix notation, are as 



follows: 



! o 




cos \p 


-sin \p 


0 




h 




^0 


= 


sin \p 


cos \p 


0 




h 








0 


0 


1 












"10 = 


ii 










T 1 




cos e 


0 


sin 6 


-J 


T 2 


h 


= 


0 


1 




0 






h 




- sin 0 


0 


cose 







and 

"21 = ' e h = 



l 2 




1 


0 


0 


h 


= 


0 


COS 0 


- sin 0 


*2 




0 


sin 0 


COS 0 



"p2 = 0i 2 = 0i p 



(3-10) 

(3-11) 



(3-12) 

(3-13) 



(3-14) 

(3-15) 



It is possible to combine the three matrices and simplify them by 
making the standard small angle approximations : 



*0 




COS 0/ 


-sin \p 


o" 




cos e 


0 


sin e 




^0 




sin 0/ 


cos xp 


0 




0 


1 


0 




\ 




0 


0 


1 




-sin 0 


0 


cos e 





0 

COS 0 
sin 0 



0 




i 






P 


sin 0 




T P 


COS 0 




k 

P 
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i 




1 


0 






io 


= 




£ 




— 0 


0 







-ip 9 

1 -0 
0 1 





*" T* ^ 




P 




T P 




k 




P 



(3-16) 



The angular velocity of the platform coordinate system with 
respect to the orientation coordinate system may be defined as: 

(3-17) 



( a ) =pi +qi + rk 

po p 4J p p 



From Eqs. (3-10), (3-11), and (3-12) we also have 



w po = w p2 + W 21 + “’lO 



= 0 ip + 0 jg + (3-18) 

In order to express Eq. (3-18) in platform coordinates, substitute 
relation for k^ from Eq. (3-12) and then make use of Eq. (3-14) as 



follows : 



u> = $ in + $ h + i£(-sin 0 i 9 + cos 0 ko) 



po 



= ? ip + [- $ sin 0 0$ cos 0 ] 



$ ip + [- i^sin 0 0 $cos 0] 






1 0 0 




i 






JP 


0 cos 0 -sin 0 




> 


0 sin 0 cos 0 




L\| 



-$sin 0 ) ip + (0 cos 0 + $ cos 0 sin 0) j p 
+ (-0 sin 0 + $ cos 0 cos 0) k 



Assuming small angles, a comparison with Eq. (3-17) gives 
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(3-19) 



p = 0 - e ^ 

♦ • 

q = 9 + 0 0/ 

r - $ - 0 § 



The rates of change of the Eulerian angles are 

• * . • 

0 = p + 0 0/ = p + 0(r + 00) « p + r0 

• * # 

9 = q - 0 0/ = q-0(r + 0 0)~q-r0 

0/ = r + 0 £ = r + 0(q- 00^) » r + 0 q 



(3-20) 



Integrating Eq. (3-20) to obtain the expressions for the Euler angles 
we obtain: 



0(t) 


= / pdt + 


f r 9 dt 




0 


0 




t 


t 


e(t) 


= / q dt - 


f r0 dt 




0 


0 




t 


t 


m 


= f r dt + 


/ q <t> dt 



o o 



(3-21) 



In the above equations, 0 , 0 , and 0/ are angles in space and when all 
are zero, the platform coordinate system coincides with the orientation 
coordinate system by definition. Having the space angles zero is no 
guarantee that the integrals of p, q and r are zero. Depending upon the 
past history of motion, p, q and r may or may not be zero. For ex- 
ample, having 0 (t) = 0 does not ensure that f 0rdt is zero. A 

o 

typical platform stabilization system can maintain the integrals of p, 
q and r at zero. Therefore, the error of the platform coordinate sys- 
tem is given from Eq. (3-21) by 
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(3-22) 



t 

A 0 (t) = f e r dt 

o 

t 

A 9 (t) = - f 0 r dt 

o 

t 

A ^/(t) = / 0q dt 

o 



If the vibratory motions are small, the integrals above become 
significant only if there is a rectification of the motions. Let 
6 = 6 Q sin cut 

r = r Q sin (cut - y ) 

then 



*0(t) = 0 o r o 



6> r 
o o 



*o r o 



/ sin cut (sin cut cos y - cos cut sin y ) dt 
o 

t 2 

f (sin cut cos y-sin cut cos cut sin y) dt 
o 

- J_ sin 2cut |cos y + -1(1 + cos 2cut) sin y 
\ 2 4 / 4 

(3-23) 



An examination of Eq. (3-23) reveals one linear term. Therefore, the 
angle drifts from null with a constant drift rate. 



A 0(t) = 



0 r cos y 
o o 

2~ 



For the system just analyzed, there will exist a coning error that 
grows with time. One possible method of reducing the error would be by 
use of a coning computer. v ' The computer would solve Eqs. (3-20) and 
compensation for the coning drift could then be built into the system. 
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3. 5 Sculling Effect 



Sculling produces an acceleration error in accelerometers 
mounted in gyro stabilized platforms undergoing combined linear and 
angular motion. This effect will be considered because what is really 
important is the effect of this error on position indication. 



The same coordinate systems and coordinate transformation 
matrix previously used in describing platform coning apply. Define 
the total acceleration vector in both coordinate systems as 



n = n xV Vp +n z k p 



n x n i o + n y j o +n z k o 



(3-24) 

(3-25) 



If there is no acceleration error, then 




n - n =0 (3-26) 

z z o 

A check on the above relations will now be made. From transformation 
(3-16) we obtain 

n = n Xo [I p - *T p + 6 V + ^ + 7 p - <t> k p J 

+ n z [ - 0 T p + ^T p + k p ] 0-26) 

Comparison of Eqs. (3-24) and (3-26) produces the relations 
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”* ■ \ + ''"Jo ‘ " \ 
n » ' 



n = » » - ■) n , r , 2 

o J o o 



and by rearranging 



n x - V = * V - 6 V 

o J o o 

n v - n v = <M Z - 
y y o z o x o 



n z ’ n z = 6 n x -<t> n y 
O 0 y o 



(3-27) 



Equations (3-27), when compared to Eqs. (3-26), show that errors be- 
tween the accelerations actually measured and the accelerations we 
desire to measure exist. 



Assume a component of the vibration of the platform is at a 

certain frequency, co. Then consider any term of Eq. (3-27), say 

0 n z . Let the linear and angular vibration vary sinusoidally at 
o 

frequency cj. 

If 0 = | 0 | sin cut 

n z = | n z I sin ( wt - y ) 




[( 1 - cos 2ojt) cos y - sin 2wt sin y ] 
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which gives a steady state drift of 




2 (3-28) 



As a consequence of the analysis, it is seen that small angular 

motions of the platform can result in rectified acceleration drifts whose 

magnitudes depend on the phase relationships between the linear and 

(91 

angular motions. The name "sculling” v 7 is given to this effect because 
the motion causing the drift is identical to that produced in propelling a 
boat with an oar at the stern. 
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CHAPTER 4 



MAGNITUDE OF THE ANGULAR VIBRATION INACCURACIES 



4. 1 General 

The errors producing forced precession of a single-degree-of- 
freedom gyroscope have been analyzed in the previous chapter. This 
chapter deals with a more practical consideration; the magnitudes of 
the various error terms. When practical ? oscillations capable of being 
produced by the Angular Vibrator , built by the Instrumentation Labora- 
tory of the Massachusetts Institute of Technology, will be assumed. 

The following are assumed as typical values for a single-degree- 

o 

- 650 gm cm 
2 

10 6 gm cm 
sec 



A(gim) — -01^45 radians 

A comparison of the various inaccuracies previously considered 
will now be made. 



of -freedom integrating gyroscope: 

AI = J IA ~ A = 
H ~ Hg P = 2 x 



C 
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4. 2 Magnitude Determination 



For convenience, the inaccuracies developed in the previous 
chapter are again listed. 

M D = (I 1A _I SA )W (IC)IA W (IC)SA + (I IA~ I SA )A (gim) (W flC)IA* W (IC)SA ) 

(3-8) 

M ge = - H SP ( w (ic)IA - A (gim) W (IC)SA ) (3 ‘ 9) 



A0( t) = / 0rdt 
o 

t 

A 6 (t) = - f 0 r dt 



A &(t) = / 0 q dt 

o 



(3-22) 



n x - V = * V - 9 V 

o J o o 

n y " n y = * V - 
J J o o o 



n z~ n z n = 6 V“^V 
o o J o 



(3-27) 



Equations (3-8) and (3-9) apply to a single gyroscope instrument while 
Eqs. (3-22) and (3-27) apply to a platform stabilized system. The mag- 
nitudes of the various terms will be determined for the case of harmonic 
angular velocities. 



4. 2. 1 Magnitude of the Dynamic Torque 



W (IC)IA = |W IA |sin(a» 1 t + 0j) 
W (IC)SA = | W SA| sin(a> 2 t + 0 2 ) 
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The anisoinertia component of Eq. (3-8), when using Eq. (4-1), 
may be written 



A I W- 



ia| 



W Sa[ [ cos K w i “ + (01-02^3 

- cos [ (o> 1 + cu 2 )t + 0 1 + 0 2 ]TJ 



(4-2) 



If 



cjf = ct> 2 = cu 



then the above reduces to: 



M 

2 





[cos ( 0 ^ — 0 2 ) ~ cos ( 2cot +01 + 02 H 

(4-3) 



Equation (4-3) consists of a constant term, and a cosine term that 
has a zero average value. Since (0 ^ - 0 2 ) represents the phase dif- 
ference between the two velocities it is evident that only the in phase 
components of the motions create an error. If the two angular veloci- 
ties are in quadrature, there will be no resultant error. 



If ^ co 2 , then the torque will oscillate at a beat frequency 
(cuj + o) 2 ) , (o>i - c Oy) and have an average value of zero. 

The second component of Eq. (3-8) for cu = may be 

written as: 



A <gim) t KaI 2 -KaI 2 + cos 2wt < I W SaI 2 cos 2 ^2 

+ |W IA | 2 cos 20 1 ) - sin2wt(|W gA l 2 sin 20 2 
+ |W IA I 2 sin 20 p] (4-4) 
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Before the magnitude can be determined, the variation of gimbal 
angle with time must be known. For the present, an ideal integrating 
gyro will be assumed so that 

A (gim) = W (IC)IA 

and 

A (gim) = / W (IC)IA dt 

-| w iaI 

= cos (wt + <f> 1 ) ( 4 _ 5 ) 

Substituting (4-5) into (4-4) we get 

|wj sin ( wt + 0 4 ) F |w lA | 2 - |w s j 2 
2 co L_ 

+ cos 2cot(|w gA l 2 cos 20 2 + lw IA l 2 cos 2 0j) 

- sin 2 cot ( (w gA | 2 sin 20 2 + (w^ 2 sin 20 (4-6) 

All the terms in Eq. (4-6) have zero as an average value; therefore, 
the entire expression may be neglected. The same applies for the 
case of ^ co 2 . 

4. 2. 2 Magnitude of the Gyroscopic Torque 

We next turn our attention to Eq. (3-9). After making substi- 
tutions for the motions we have 

M ge = ~ H SP t |W li sin ( Wl t + 0 1 ) - A (gim) |wj sin (co 2 t + 0 2 ) ] 

(4-7) 

The first term in the bracket will have a zero average value. Our 
attention will now focus on the cross coupling coning term: 



32 



A (gim) H SP KJ sin < ‘V + <j> 2 ) (4-8) 

Again assuming a perfect gyro so that relation (4-5) applies we obtain 

- H spl wj | W gA l cosCwjt+^j) sin(cu 2 t + 0 2 ) (4-9) 

Cl) 

o) 1 = cug = ci) 

we can write: 

' ~ H SP I W Ia| |. W sa| [ sitl (0 2 “ 0 l) + sin ( 2wt + 0 2 + 0 1^ 

2 0 ) 

(4-10) 

Thus the torque depends on the out of phase components of the veloci- 
ties. Should the two velocities be in phase there would be no torque. 
Only velocities in quadrature produce a coning error. From this 
discussion it is seen that the float drifts from null at a constant rate 
of 

-t-KJ |w d sin(02“0^) rad/sec (4-11) 

2cu 

For the condition of cu^ ^ Eq. (4-9) contains no linear or 
average torque values. 

Note that when coning is a maximum, the anisoinertia effect is 
a minimum. Angular motions that are not exactly in or out of phase 
will produce both effects. However, even for this condition it is seen 
that the anisoinertia error may be neglected, since it is very small 
in comparison to the cross coupling coning error. This is due to the 
fact that the coning term involves Hgp which is very much larger than 
any of the values in Eq. (4-3). 

The coning drift rate is plotted in Fig. 4-1 as a function of 
frequency for angular velocities 30, 60 and 90 degrees out of phase. 
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DRIFT RATC 



All of the previous work has been for an ideal gyroscope. Ac- 
tually, the true gimbal angle variation should be obtained from the gyro 
equation 



A (gim) + W (IC)SA A (gim) " W (IC)IA 



(4-12) 



If Eq. (4-12) is used to obtain the gimbal angle motion, the only differ- 
ence from the previous result would be a nonzero average value of coning 
drift. This can be easily shown. First integrate Eq. (4-12) 



A (glm) ■ * J ^ /; w (IC)IA « 

(4-13) 

Assuming sinusoidal variations and small magnitudes for the angular 
velocities, the exponentials may be approximated by the first two terms 
of a series 

± f W (IC)SA dt i . r w Ht 

“ 1 ± / W (IC)SA dt 

(4-14) 



Using the 



mg me relations (4-1) and the gimbal angle expression 
can determine the nature of the coning term 

A (gim) " / IV ! W Sa' — sin(« 1 t + ^ j ) cos(w, 

o Ct>2 



o 

-J- iw J 1 w sa ! j 
2 



(4-13), we 



^t + 0 2 ^ 



+ <j >l - w 2 t -<j) 2 ] 



sin [oijt + Wgt + <p i + 0 2 ]] 



dt 



if = w 2 



0J o = CJ 



then 

A (gim) ~ “ > W IA' I W Sa' sin ^ ^ 1 _ ^ 2' 



cos(0 1 +<j> 2 , 

,) xt + i 

2 (jJ 



y) - cos(2o>t+ 
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It is seen that a constant drift rate appears similar to expression 
(4-11) plus an oscillatory term and a nonzero average term. Since 
the constant drift rate term dominates, the approximation for gimbal 
angle motion using Eq. (4-5) is valid. 

4. 2. 3 Magnitude of Platform Coning 

From Chapter 3 it was seen that platform coning effect causes 
the Euler angles to drift from null at a rate equivalent to 0(t) 

0 I* 

A?(t) = 0 0 cos y (4-15) 

2 

The above equation was the result of integrating one of the error equa- 
tions (3-22). 



If it is assumed that 

r = r Q sin (cot + y ) « $ 
i// = i// Q cos(cot + y) 



then 



^ = 



and Eq. (4-15) is now 



- co sin (cot + y ) 

-V 



A?(t) = 



- '/'o 8 o w cos y 
2 



(4-16) 



The assumption will now be made that the amplitudes of the 
motions are known but that the directions and phase relationships be- 
tween the components of the motion are random. 

Let 

'p Q = K cos /3 
0 Q = K sin (3 
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Hence 



a$(t) = 



sin 2 /3 cos 7 



K 2 a) 



Now assume that /3 and y are independent and uniformly distributed 
from 0 to 2tt, so that the mean value of each of the above trigonometric 

terms is zero, and the variance of each is 1/2. For example 

2ir 

Mean (cos y) = (jl cos y = E(cos y) = f cos y dy = 0 

2tt o 



2 2 i 

Variance (cos y) = cr (cos y) = E (cos 7) = 

2ir 



r 2 * 2 

J cos y dy - 
o 



1/2 



We can now write 



2 

2 

E [ sin2j3 cos 7 ] 





and finally the standard deviation or RMS value is 



a 



K 



cu 



(4-17) 



where 

K = maximum amplitude, radians 
c 0 = frequency, rad/sec 

Equation (4-17) is plotted in Fig. 4-2. 

4. 2. 4 Magnitude of Sculling Effect 

The sculling effect will be treated by the method used in the 
previous section. Consider the first of Eqs. (3-27) : 

n „ - n „ = f n - 9 n (4-18) 
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Using the form of the sculling error found in Chapter 3 we can write 



n x" 


\ ■ T 1*0 \J “ 5y i -? 


K n zJ cos 


If we let 






K.I 


= n cos 0 




'"‘o' 


= n sin /3 




1 *0 1 


= K cos a 




1 e o 1 


= K sin a 




then 






n x - n x 


- — nK (cos a cos /3 cos y 1 

n 9 1 


- sin a sin /3 



(4-19) 



(4-20) 



Consider that a, 0, y^, and y 2 are independent of each other and 
uniformly distributed from 0 to 2 ir. The mean value of each trigono- 
metric function is zero and the variance of each is 1/2. Then 

E [(n x ~ n x ) 2 ] = (— nK) 2 E[ cos 2 a cos 2 /3 cos 2 y 1 + sin 2 a sin 2 /3 cos 2 y 2 ] 
o 2 

1/2 

~ n x > = - nK [ 2 ( X) (1) (i.) ] 
x x o 2 2 2 2 

= standard deviation or RMS (4-20) 

4 



where 2 

n = maximum acceleration (ft/sec ) 

K = maximum amplitude (radians) 



Equation (4-20) is plotted in Fig. 4-3. 
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NOlJ-OW jianj-ntu/i/ 



CHAPTER 5 



THE RESPONSE OF THE SINGLE -DEGREE -OF-FREEDOM 
INTEGRATING GYROSCOPE TO TWO ANGULAR VIBRATIONS 
SIMULTANEOUSLY APPLIED ABOUT THE SPIN REFERENCE 
AND OUTPUT AXES 



5. 1 Introduction 

Gyroscopic effects caused by sinusoidal angular vibrations are studied 
in this chapter. Two vibrations of different frequencies have been applied to 
two of the body fixed axes of the gyroscope case. The particular axes to which 
these inputs are applied are the output and the spin reference axis. 

A series solution has been used to approximate the theoretical angle which 
should exist between the float and inertial space about the output axis. The 
method used was applied to a case of two equal input frequencies in ref. 10. 
Thus, this is the more general situation and a further study of the problem with 
varying inputs. 

A typical dynamic model of a gyroscope has been used as the basis for 
the equation of motion. Effects such as inertia damping, and compliance have 
been included. The resulting equation presently has no closed form solution 
and the approximate series method was used. 

The investigation was an attempt to discover if and when gyroscopic drift 
can be expected along with the associated values of magnitude. . 
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5. 2 Definition of Symbols 



The following symbols and definitions have been used in the analyses 
(units are in parentheses) 



a 

b 

c 

c d(OA) 

h 



k 

t 

x 

A (I-C)NS 



a (i-f)oa 

A (I-C)OA 

A 

B 

C 

D 

E 

F 

G 



= constant (radian) 

= constant (radian) 

= constant (radian) 

= damping coefficient between float and case about the 
output axis (dyne -cm -sec/radian) 

H 

c d(OA) 

K (C-F)OA 

c d(OA) 

= time (seconds) 

= the variable representing 

= angular rotation of the case about a stationary axis 
aligned with the earth's axis which coincides with the 
spin reference axis when A^ qjqA = 0 (radian). 

= angular rotation of the float about the output axis (radian) 
= angular rotation of the case about the output axis (radian) 

= constant (radian) 

= constant (radian) 

= constant (radian) 

= constant (radian) 

= constant (radian) 

= constant (radian) 

= constant (radian) 
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H 



I 



H 1 

(F)OA 
I 



K 



W 
J 
K 

(C-F)OA 



angular momentum of the wheel (dyne-cm-sec) 
constant (radian) 

inertia of the float about the output axis (gm-cm^) 
inertia of the wheel about a radial axis (gm-cm^) 

constant (radian) 
constant (radian) 

elastic restraint about the output axis between float and 
case. Rotational stiffness between float and case about 
the output axis (dyne-cm/rad) 



L 

M 

N 

"1 



constant (radian) 
constant (radian) 
constant (radian) 

circular frequency of vibrational input (rad/sec) 
circular frequency of vibrational input (rad/sec) 



^EQOA + *W 
c d(OA) 



gyroscope time constant (seconds) 



(The inertia of the wheel is added to that of the float in r 
since the wheel is assumed to be perfectly rigid with 
respect to the float about the output axis in the analyses.) 



5. 3 The Equation of Motion of the Dynamic Model^ 

The simplified equation of motion of the dynamic model shown in Fig. 5-1 



is given by 



K 



(C-F)OA 



A (I-F)OA - A (I-C)OA 



] + c d(OA) [ A (I- 



F)OA * A (I-C)OA 



■] 



! (F)OA + 1 



wj A (I-F)OA ■ H J^ A (I-C)NS A (I-F)OaJ (5-1) 
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By using small angle theory, the term H |^ A (i_c)NS A (I-F)OaJ is introduced 
by the component of angular velocity about the input axis produced by the motion 
about the space-fixed axis. Division by and rearrangement yields: 



K (C-F)OA 


• 

A , A , 


Vjoa + *w 


c d(OA) 


(I-F)OA + "(I-F)OA + 


c d(OA) 



(I-F)OA 



+ 



-H 

_ c d(OA) 



A (I-C)NS A (I-F)OA 



K (C-F)OA 

c d(OA) 



A (I-C)OA + 



A (I-C)OA 



or 

j^k - hAg_ C ) NS J A (x_y)0A + A (I-F)OA + 

= A (I-C)OA + 

The stiffnesses between the wheel and float about the input and output axes are 
assumed to be infinite compared to the inertia of the wheel for the frequency 
range considered. The stiffness and/or damping between the float and case 
about the input axis is also assumed infinite compared to the inertia of the 
float and wheel within the frequency range considered. 



7 A (I-F)OA 

kA (I-C)OA ( 5 ' la) 



5. 4 Method of Solution of the Equations of Motion for Two Different 
Sinusoidal Vibration Inputs of Small Magnitude 

The inputs introduced to the gyroscope case may be represented as follows: 10 



(I-C)NS = a sin “l 1 


(5-2) 


(I-C)OA = b sin w 2 t + c cos w 2 t 


(5-3) 



In most applications of the single-degree-of-freedom integrating gyroscope, 
the constant h is designed equal to unity. For purposes of simplification, this 
constant will be included in the vibration inputs as a factor of the constant a. 
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the SDFI Gvroscope for Motions 
about the Spin and Outwit Axes. 



To be strictly correct, earth rate should be included in the definition of the 
vibration input about the North-South axis. This rate is, however, small 
compared to the input angular velocity. 

Substitution of equations (5-2) and (5-3) into equation (5-la) yields a 
linear differential equation with a variable coefficient of the following form: 

-acu^ cos cu^tj x + x+ rx 

= bc ^2 cos cugt - cu >2 sin cu 2 t + kb sin o> 2 t + kc cos a> 2 t (5-4) 

where A (IF)0A = x 

There is no known closed form solution for equation (5-4). However, if the 
constants are small compared to unity, a suitable series solution may be ob- 
tained. * 

5.4.1 Negligible Elastic Restraint Between the Float and Case 
About the Output Axis 

When the elastic restraint is very small compared to the damping con- 
stant, i. e. , 

k -> 0 

equation (5-4) reduces to: 

(-acuj cos cu^t) x + x+ rx = bo> 2 cos o> 2 t - ccu 2 sin cu 2 t (5-5) 

Assume a solution of the form: 

x = Xj + x 2 (5-6) 

where x^ is defined by: 

Xf + TXj = bo> 2 cos co 2 t - cousin u> 2 t (5-7) 

Substitute (5-6) and (5-7) into (5-5): 

(-acoj cos cu^t)(x 2 + x 2 ) + (x^ + x 2 ) + r(x^ + x^ = bu> 2 coscu 2 t - co> 2 sincu 2 t 

(5-8) 

♦The series solution technique was suggested by Mr. Daniel Goldenberg 
of the Mathematics Group, MIT Instrumentation Laboratory. 



or 



(-acu^ cos Wjt)x2 + Xg + tx ^ = (acu^ cos cj^t) 
Now let x^ = x^ + x^ 
where is defined by: 

x 3 + rx' 3 = (ao> 1 cos cu^t) x^ 
Substitute (5-9) and (5-10) into (5 -8a): 



(5 -8a) 
(5-9) 



(5-10) 



(-acu 1 cos cu^tKxg + x^) + (x^ + x^) + t(x^ + x 4 ) = (acuj cos co^t) x ^ (5-11) 



or 

(-aa>i cos cu^t) x^ + x^ + rx^ = (acu^ cos cu^t) x^ (5-lla) 

Continuation of this process results in a series solution of the following form: 

x = xi + x 3+ x 5 + ••• x 2m-l + x 2m (5 ' 12 > 

where X 2 m is the even term remaining after the above process is repeated 
m timeSo X 2 m satisfies the equation: 

(-a^i cos w it) x 2m + x 2m + TX 2m = (aw l cos x 2m-l (5 - 13) 



Problem: 

(-acu 

let 

where: 







A Brief Summary for 


k 0 is: 




Wjt) : 


X + X 


+ TX = bcUgCOSCUgt- 


ccug sin cjgl 


(5-5) 


X = 


x l + 


x 3 + x 5 + . . . . x 2m 


-1 + X 2m 


(5-12) 


x l + 


•• 

TXj 


= be u 2 cos cu 2 t - co >2 sin o> 2 t 


(5-7) 


x 3 + 


•• 

rx 3 


= (aa>i cos cu^t) x^ 




(5-10) 


x 5 + 


tx 5 


= (aaj^ cosa^t) x 3 




(A -6) 


• 

X (2N+1) + 


tx ( 2N+1) = ^ os ^ 


1^ X (2N+1) 


(A -6a) 
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Each successive term in the series of equation (5-12) can be shown to be 
smaller than the preceding term by a factor of approximately a. (Refer to 

o 

ref. 10 for further discussion.) Since a is of the order of 2 x 10” radians, 
this series is seen to converge rapidly. Hence all terms beyond x^ can be 
reasonably neglected. 



The following evaluation of the series and associated coefficients clearly 
shows the appearance of the factor a which permits terms beyond a certain 
point to be neglected. (See Appendix A for solution details.) 



A (I-F)OA “ x 



k = elastic restraint = 0 



x = x i + x 3 + x 5 + °ther terms which are neglected 



x^ = A sin o> 2 t + B cos u> 2 t 



Xg = C sin cogt + D sin cu^t + E cos cu^t + F cos cu^t 



(U 3 = CUi + cu 2 



= - w 2 



(5-12) 
(A -7) 
(A-17) 



x 5 = G sin (2cu^t + cu 2 t) + H 1 cos(2a> 1 t + a> 2 t) + J sin (2^^ - a> 2 ) 

+ Kcos (2u>jt - a> 2 t) + Lsin(-cu 2 t) + Mcos(-a> 2 t) 

+ Nsin(cu t) + Pcos(cu 9 t) (A-37) 

2 z 



A = 



b + c to> 2 

~ 2 2 
1 + T CU 0 



(A- 12) 



act?! (B - t Acu q) 

C = i (A -24) 

2cu 3 (1 + r 2 co 2 ) 



-acu^A + rBcu Q ) 

E = i — (A-26) 



c- bro> 0 

B = tL (A- 14) 

i 2 2 

1 + T Ct> 2 



acu^B + tAcu.) 
D = 1 i- 

2cu 4 (1 + t 2 cu 4 ) 

acu^A - tBcu 4 ) 
F = 1 q 



(A-28) 



(A -30) 



o/i 2 2. 

2cu 3 (1 + r a> 3 ) 
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o/i 2 2i 

2w 4 (1 + r w 4 ) 



— 3.(^2 [C + T E(2(x)^ + CO 2 ) ] 



G = 



ao> 4 [E - rC^cu^ + o> 2 ) ] 
2(2oj^ + o> 2 ) [1 + t 2 (2o>j + o > 2 ) 2 ] 



(A- 47) 



H l = 



2(2w 1 + w 2 )[l + t 2 (2 cj 1 + w 2 ) 2 ] 



(A- 48) 



act) 1 [F - t D( 2o? 1 -o) 2 ) ] 
2(2ct)^ - o> 2 ) [1 + t 2 (2o^ - o> 2 ) 2 ] 

(A -50) 



aa)^ [- E + tCo> 2 ] 
2ct> 2 [1 + r 2 a? 2 ] 



(A-51) 



ac^ 1 [F + r Do> 2 ] 
2oj 2 [1 + t 2 cu 2 ] 



(A- 53) 



-acoj [D + rF(2a> 1 - ct> 2 ) ] 
2(2w 1 - o> 2 ) [1 + r 2 (2a) 1 - c 2 ) 2 ] 

(A -49) 



-aco 1 [C + tEcl ) 0 ] 
M = fL 

2 0> 2 [1 + t 2 Ct) 2 ] 



(A- 52) 



aw 1 [D -rFw 0 ] 

P = — i i— (A- 54) 

2o> 2 [1 + r 2 a) 2 ] 



r 



max 
:-Jb 2 +c 2 j 



max 



. 0009 sec = . 0 X 10 ^sec 
2 x 10“^ radians 

4 X 10" ^ radians 



> 



Typical values as used by 
the MIT Instrumentation 
Laboratory Angular Shaker 



With negligible elastic restraint between the float and case, the series 
solution consists of sinusoids. Such sinusoids will exhibit a time-average 
value of zero except for special cases which are of interest. 

The first case has been studied in ref. 10 and concerns equal frequencies, 

W 1 = w 2 

As in further explained in Appendix A, equal input frequencies will cause the 
terms of D sin o> 4 t and F cos cu^t to have indeterminate values. By applying 



L' Hospital's Rule, a constant drift term which is time dependent is found 
to exist. This term is: 



Drift = t where cv = = cog 



[■ 



where co = = cog 



This existing drift is identical to that found in ref. 10. 

The second case of interest concerns the almost natural region for sus- 
picion, namely harmonics. Terms concerning x^ of the series give indetermi 
nate values when 2co ^ = u> 2 ° By L' Hospital's Rule, a drift term appears which 
is very similar to the term arising from equal frequency input. This term is: 



If further series terms such as x^ and Xg were evaluated, drift terms 
associated with other harmonic multiples are found. Thus, the conclusion is 
that equal frequencies and harmonics are capable of causing drift for the case 
of negligible elastic restraint. 

The physical case will probably display many superimposed input fre- 
quencies. This is likely to have some harmonics thereby creating drift. Such 
drift is substantiated by this theoretical solution but it is suspected that fre- 
quencies approaching the harmonics may well create actual physical drift. 

5. 4. 2 Finite Elastic Restraint Between the Float and Case 

About the Output Axis 

For the case when k 4 0 , it is desired to find the solution of the 
following equation: 

(k - awj cos cu^t) x + x + rx = ba >2 cos - cc*^ sin Wgt + kbsincugt + kccoscu 2 t 




See Appendix A for proof and 
necessary mathematics 



(5-4) 
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The mathematical solution with details is in Appendix A. The resulting series 
with values of coefficients in terms of the known constants is as follows: 



A (I-F)OA ■ x 



k / 0 



x = + Xg + x 5 + other terms which are neglected 

x 1 = A sin cu 2 t + B cos a^t 



(5-12) 

(A-67) 



Xg = D sin cu^t sin u^t + E sin cu^t cos cu 2 t 



+ F cos co^t sin cu 2 t + G cos cu^t cos u> 2 t 



(A- 7 5) 



Xg = sin 2cu^t sin cu 2 t + I sin 2cu^t cos cu 2 t 
+ J cos 2a>it sin cu 2 t + L cos 2cu^t cos a> 2 t 
+ M sin u> 2 t + N cos cu 2 t 



(A-89) 



Expressions for the coefficients in terms of known constants may be found in 
Appendix A. The coefficients have been left in determinant form. 

By the use of suitable trigonometry relations the following forms are 
obtained for x^: 



When = a> 2 , the case of ref. 10 occurs. An indeterminate form does not 
appear because k is finite. Thus time dependent drifting is not a factor with 
non- zero elastic restraint. However, this special case of equal frequencies 
gives a value of unity to terms containing cos(cu^ - cu 2 )t and the result is a 
steady non -zero value for the angle between the float and inertial space about 
the output axis. This is verified by ref. 10. 




51 



The other case of interest with finite elastic restraint again concerns 
harmonics. This can easily be seen by converting terms by means of 
trigonometry relations to the following: 



x 



5 



+ 



+ 



+ 




2 



2 

J_ 

2 

2 



[cos(2ct? 1 - a> 2 )t] 

[sin^o^ + w 2 )tl 
[sin(2w 1 + o> 2 )t] 

[cOS^CUj + w 2 )t] 



^1 [cos^a^ cu 2 )t] 
2 

+ -L [sin(2co 1 - cu^t] 

- -y [sin(2co 1 - w 2 )t] 

+ — [cos(2co 1 - ^ 2 )t] 
2 



+ M sin a> 2 t + N cos cu 2 t (A-89a) 

The drift resulting from the harmonic when - a> 2 is clearly seen as 
a possibility. However, because of a finite value for elastic restraint, indeter- 
minate forms do not appear. Nevertheless, this particular condition causes an 
angle to take place about the output axis between the float and inertial space. 
This occurs when cos(2cOj - a> 2 )t is unity in value. 



This case gives a term showing that a constant angular displacement 
occurs for the harmonic condition instead of a time dependent drift term. 



Such a condition is almost a natural situation for a gyroscope. As the 
finite restraint is increased, the displacement between the float and the refer- 
ence in inertial space should become smaller. The extreme limits are for 
zero elastic restraint and infinite elastic restraint. The former case gives a 
time dependent drift while k = oo creates a situation permitting no angle, what- 
soever, to exist between the float and inertial reference. 

The most interest undoubtedly lies with the case of finite, non- zero, 
elastic restraint. This is the existing case of practical importance. The 
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concluding results are that 1) a time -independent drift term is possible if 
input frequencies contain harmonics or equal values, 2) float angular dis- 
placement with reference to inertial space depends on the degree of elastic 
restraint in the gyroscope, and 3) high harmonic terms may possibly be 
neglected because of their small value which becomes insignificant when 
compared to manufacturing tolerances. 

It is pointed out that a drift may appear to be present at any one instant 
of time if a time -averaged measurement is not obtained even if harmonics are 
not present. This could be of noticeable magnitude dependent upon the input 
frequencies and their own absolute magnitudes. 

It is strongly recommended that this theory be applied to the MIT Instru- 
mentation Laboratory angular shaker in an attempt to verify the results. Har- 
monic inputs appear feasibile for the shaker by minor adaptation after sufficient 
data has been collected concerning the case of two inputs of equal frequency. 

Terms in the approximate series solution beyond have only theoretical 
interest at present. Until gyroscopic tolerances exceed their present values, 
common sense tells us that the selected series is as good for analysis as would 
be the unknown closed form solution for the equation of motion. 

5. 5 Drift of a Typical Gyroscope for the Case of Negligible Elastic 

Restraint About the Output Axis with Two Equal Input Frequencies 

The ensuing results are based on the series solution for the equation of 
motion when the elastic restraint between the case and float about the output 
axis is negligible. A further stipulation is that = cog in order that this in- 
formation may apply directly to the MIT Instrumentation Laboratory angular 
shaker. 
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The angular shaker has a frequency range from twenty (20) cycles per 
second to one hundred (100) cycles per second. The constant "a” has a max- 
imum value of . 002 radians and the constants ”b” and ”c” are limited such 

I 2 2 * 

that ^ b + c has a maximum of ,004 radians. A typical value of r, the 
time constant of the damper-mass system, is .0009 seconds. 

By definition: 

A (I-C)NS = asinw it (5-2) 

A(i_c)oa = b sin CJ 2 t + C cos a) 2 t (5-3) 

For this particular case, cu = cu^ = as the angular vibrator is capable of 
applying only one frequency at a time to the two input axes. 



If c = 0, a (j_c)ns and A (l_C)OA are in Phase and they will differ only 
by the magnitudes of "a” and ”b”. (a = b is a special case.) In general, "b" 
and "c” will have finite values thus causing a (i_q)nS anc * ^(I-C)OA to 
by a phase angle as well as magnitude. (Here, again, a = tF + c^ is a 
special case.) 

b sin ( ji ) t + c cos cut = sin(cut + cot"^b/c) 






. 2 2 
b + c 



(5-14) 



cot"*b/c 



= 0 , a phase angle 



If cu^ = cog, the following is drift rate for k 0 : 



Drift Rate = 

2 



1 T CU 

C 

i 2 2 

1 + T CU 



A plot is made for drift rate/ ac for various ratios of b/c. 



Drift Rate 
ac 



= jrf 



1- — rcu 
c 



- 2 2 

1+7 CU 



1 T CJ 

C 

- 2 2 
1 + T CU 



(5-15) 



(5-16) 
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Positive and negative values of b/c were selected to correspond to phase angles 
from approximately ± 5° to ± 84°. Ratios for b/c selected are indicated on the 
graph. 

The curves of Fig. 5-2 show that for a frequency range of 20 to 100 cycles 
per second, the drift rate will undergo a change in direction for b/c = + 2. 0 
around 90 cycles per second. The frequency required for this reversal is fur- 
ther decreased for increasing values of b/c. 

Although the curves represent specific values for constants, it is believed 
to be typical of gyroscopes which may be analyzed by the MIT angular shaker. 

It is suggested that these curves be compared to actual data from the vibrator 
test machine at the earliest opportunity. 
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drift rate radios 

Ouc m SECOND 



CHAPTER 6 



CONCLUSIONS AND RECOMMENDATIONS 



6. 1 Conclusions 

1. When a floated single -degree-of-freedom integrating gyro- 
scope is subjected to vibrations undesired torques or drift 
rates may exist. A summary of such unwanted errors, 
arising when the vibrations are angular in nature, is pre- 
sented in table form. The errors are classified by name, 
type and cause. 



Name 



Type Error 



Anisoinertia effect Dynamic 



Cross coupling coning Gyroscopic 



Cause 

Simultaneous angular 
motions about the case 
input and spin axes 
along with unequal mo- 
ments of inertia about 
these axes. 

Float precession about 
its output axis which 
causes a component of 
the case spin axis angular 
velocity to be coupled to 
the float input axis angu- 
lar velocity. 
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Simultaneous angular 
velocities about plat- 
form axes such that 
the float input axis 
undergoes a periodic 
conical motion 

Simultaneous angular 
motion about any plat- 
form axis and linear 
acceleration along a 
perpendicular axis 
producing a rectified 
acceleration along an 
axis perpendicular to 
both axes. 

2. A second table is presented to illustrate which effects are 
present when a single -degree -of -freedom integrating gyro- 
scope is subjected to constant angular velocities and periodic 
vibrations which may be in or out of phase. 



Constant Angular 
Velocities 


Periodic In Phase 
Angular Velocities 


Periodic Out of Phs 
Angular Velocities 


Anisoinertia 


Yes 


Yes 


No 


Cross coupling 


Yes 


No 


Yes 


coning 


Platform coning 


Yes 


No 


Yes 


or gimbal walk 


Sculling* 


Yes 


Yes 


No 



Platform coning or Geometric 

gimbal walk 



Sculling effect Geometric 



♦For sculling the phase refers to the phase relation between 
the linear acceleration and the Euler angle motion (space angle motion). 
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3. Coning is a cause of drift for a single -degree -of -freedom 
integrating gyroscope, 

4. Cross coupling coning is inversely proportional to frequency 
and directly proportional to the product of the amplitudes of 
the angular motions. 

5. Cross coupling coning may be minimized by use of a tight 
servo loop which keeps the gimbal angle close to zero. 

6. Platform coning is directly proportional to frequency and di- 
rectly proportional to the square of the amplitude of angular 
motion. 

7. Sculling generated drift is directly proportional to the ampli - 
tude of the angular motion and linear acceleration . 

8. For zero elastic restraint between the float and case about 
the output axis: (Based on a series solution of a dynamic model) 

(a) Constant drift rate appears when two equal frequencies 
are simultaneously applied to the spin and output axes 
(e.g. , Wj = « 2 ) 

(b) Constant drift rate appears when the frequencies applied 
simultaneously to the spin and output axes have a har- 
monic relation (e.g. , 2cjj = 

(c) The total drift at any instant of time consists of several 
superimposed sinusoids. Except for the special cases 
concerning inputs of equal or harmonic frequencies re- 
ferred to in (a) and (b), the time averaged value of 
angular displacement of the float with respect to inertial 
space about the output axis is zero. 
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9. For finite elastic restraint between the float and case 
about the output axis: (Based on a series solution of a 
dynamic model) 

(a) Drift terms which are time independent appear 

when two equal frequencies are simultaneously 
applied to the spin and output axes (e. g. , = a> 2 )* 

These drift terms have a constant magnitude de- 
pendent upon a, bj c, r, (Up and but not 
dependent upon time. 

(b) Drift terms which are time independent appear 
when the frequencies applied simultaneously to 
the spin and output axes have a harmonic relation 
(e. g. , 2cj^ = (Ug)- These drift terms have a con- 
stant magnitude dependent upon a, b, c, r, o>p 
and oj 2 but not dependent upon time. 

(c) The total drift at any instant of time consists of 
several superimposed sinusoids. Except for the 
special cases concerning inputs of equal or har- 
monic frequencies referred to in (a) and (b), the 
time averaged value of angular displacement of the 
float with respect to inertial space about the output 
axis is zero. 

6. 2 Recommendations 

1. Verify theoretical results by use of the MIT precision angular 
vibrator. This should be done for various values of elastic 
restraint. 
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2. The present vibrator is capable of applying two simultaneous 
inputs of equal frequency to two orthogonal gyroscope axes. 

It is suggested that the machine be investigated to determine 
if input frequencies with harmonic relationships can be simul- 
taneously applied. 

3. Linearity should permit results obtained from a specific fre- 
quency applied to one axis to be superimposed upon the results 
obtained from a different frequency applied to another gyroscope 
axis. With th^s procedure, it may be possible to obtain physical 
data concerning simultaneous inputs of any desired frequencies 
within the limit of the vibrator. Information for assorted com - 
binations of input frequencies to two orthogonal gyroscope axes 
could thus be obtained from an angular vibrator capable of 
applying only one frequency at a time. 

4. Great care should be taken to obtain valid data in the low fre- 
quency range. This is the suspected source of many uncertainties 
and may be caused by the relatively obscure bending modes of 
large missiles and high speed aircraft. 
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APPENDIX A 



SERIES SOLUTION FOR THE EQUATION OF MOTION OF A 
SINGLE-DEGREE-OF-FREEDOM INTEGRATING GYROSCOPE 



A. 1 Series Solution for Negligible Elastic Restraint Between 
the Float and Case About the Output Axis (k->*0) 

The general equation is: 

(k - aa>i cos co^t) x + x + rx = bc^ cos - ccug sin o^t + kb sin 

+ kc cos cugt (A-l) 

Ifk->0, we have: 

(-acuj cos cu^t) x + x + rx = bcug cos c Ugt - ccug sin cUgt (A-2) 

In accordance with the series established in Section 5. 5, the interest lies with 
the following terms: 

x = xj + x 3 + x 5 + . . . . + x 2m -i + x 2m 

where: 

Xf + rxj = bu) 2 cos cu 2 t - ccug sin o> 2 t 
x 3 + rx 3 = (a coj cos cj^t) 
x 5 + tx 5 = (acuj cos cujt) x 3 

X (2N+1) + TX (2N+1) = (aw l cos w l t)x (2N+l) 

Let 

Xj = A sin Wgt + B cos Wg t 
Xj = A o >2 cos Wgt - B a >2 sin a^t 
Xj = - A sin - B a)| cos Wgt 



(A -3) 

(A-4) 
(A-5) 
(A -6) 
(A-6a) 

(A-7) 
(A- 7a) 
(A- 7b) 
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A >- A 



Substitute (A-7), (A-7a), and (A-7b) into (A-4) 

2 2 

[Aoj 2 cos Bcu 2 sin 00 2 ^ + 7 [~ Aa, 2 sin ^2* ” Bct, 2 C0S 



= bcog cos ci^t - ca >2 sin c^t (A-8) 

Equate coefficients of like terms and solve for A and B. 

o 

Aci^-TBcUg = bo >2 (from cosine terms) (A-9) 

2 

-Ba> 2 -TAa >2 = “ cw 2 (from sine terms) (A-10) 

A = b + tB o >2 (A-9a) 

B = c — t A co g (A-lOa) 

A = b + rcu 2 [c - TAcug] (A-ll) 

2 2 

A + t o^A = ^ + crCt, 2 



b + CTCJ 0 

A = — 

i 2 2 

i + T 0) 2 

B = c-rci) 2 [b+ rBo^] 

B + B = c - br o> 2 

c - br cu 0 
B = 

- 2 2 

1 + r o> 2 

Only those terms affecting the steady state solution are considered. 
An equation for is obtained by substituting (A-7) into (A- 5). 
+ tx 3 = (acuj cos (UjtMA sin cjgt + B cos o^t) 

• •• 

x 3 + TXg = Aaw^ cos Wjt sin o> 2 t + Baw^ cos o>jt cos u> 2 t 



(A- 12) 



(A-13) 



(A- 14) 



(A- 15) 
(A -15a) 
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From basic trigonometry relations, the following are obtained: 

cos cu^t sin = — [sinCu^t + a^t) - sin(co^t - a^t) ] 

2 

COS CJjt COS Cdgt = — [cOs(cUjt + C^t) + C0s(cj^t - COgt) ] 

Let cjgt = (u^t + u> 2 t) and = (o^t - cu 2 t) 

In view of these relations: 

t ## Aaw 1 Bacu-, 

x 3 + rx 3 = t sin w 3 t ” sin w 4^] + [ cos w 3 t + cos ^4^1 (A-16) 

2 2 

Let 

Xg = C sin cjgt + D sin cu^t + E cos cu^t + F cos cu^t (A- 17) 

Xg = Ccjg cos ajgt + Dcu^ cos cj^t - Eajg sin w^t - Fcl> 4 sin cu^t (A- 17a) 

•• 2 2 2 2 

x 3 = -Co> 3 sin cj^t-Do ^ sin co 4 t-Ecu 3 cos cu^t-Fcu^ cos co^t (A- 17b) 

Substitute (A- 17), (A-17a), and (A-17b) into (A-16) . 

[Ccug cos cugt + Dcu 4 cos cu^t - Ew^ sin cu^t - Fcj 4 sin cu^t] 

2 2 2 2 
+ [-tCco 3 sin cu^t -rDcu| sin o^t -rEcu^ cos cu^t -tFcu 4 cos cu^t ] 

Aacu^ Aacu* Bacu.. Bacu^. 

= sin cu Q t sin oj A t + cos cu Q t + cos wA (A-18) 

2 6 2 q 2 6 2 



Equate coefficients of like terms and solve for C, D, E, and F. 

2 Bac^ 

cos cu 3 t : Cct> 3 -rEcu^ - 



cos cu^t : Do> 4 -tFcu 4 



sincu^t: -Eco 3 -rCo> 3 



2 

Bacu^ 

2 

Aacu., 



(A- 19) 
(A -20) 



(A-21) 
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sin o)^t ; 



-Fcu 4 -tDcd 4 = 



Aacu-, 



C = 



o 

1/2 Baa^ + rEa)^ 



E = 



-1/2 Aaa^ - rC cu 2 



Substitute (A- 2 la) into (A- 19) 



Co> 0 -TCUo 



-1/2 Aacu^ -rCco^ 



Baco., 



Co> 3 + l^rAacu^cUg + t 2 Ccu 3 = 1/2 Bao^ 



C(o> 3 + t 2 cu^) = 1/2 Bacu^ - 1/2 rAao^cUg 




Substitute (A- 19a) into (A- 21) 



-Ecu 3 - rcu 3 



1/2 Bacu^ + tEcu 3 



CU Q 



Aacu., 



-Ecu^ - 1/2 rBacu^ cu 3 -t 2 Ecu 3 = 1/2 Aacu^ 



9 o 

E((l)^ + T Wg) “ — l /2 Asuij - 1/2 T B&WjWg 



E = 


- Aaa>i - r Bacu^cug 




2(o» 3 + r 2 ) 



(A-22) 
(A- 19a) 

(A-21a) 

(A-23) 

(A-24) 

(A-25) 

(A-26) 
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D 



1/2 Baa>i + tFo)^ 



(A- 20a) 



F = 



1/2 Aacj^ - rD 



< 



Substitute (A-22a) into (A-20) 



Dcu. - TO). 



1/2 Aacuj - rD u> 4 






Bacu-, 



p q 

Do> 4 - 1/2 rAac^ct^ + r Do>| = l/2Bau> 1 
D(a> 4 + r 2 cuj = 1/2 Bacu^ + 1/2 rAacu^cu^ 




Substitute (A-20a) into (A-22) 



-Fcu 4 - T 0> 4 



l/2Bao>^ + tFo> 4 



- 1/2 Aao> 1 



-Fct> 4 - 1/ 2 r Baa^u> 4 - t 2 Fo> 4 = - 1/2 Aacu^ 



F(o> 4 + r 2 o> 2 ) = 1/2 Aao^ - 1/2 rBaa^a^ 



F = 



Aao)^ - TBacu^cu 4 

2 (o> 4+ r 2 w |) 



An equation for x & is obtained by substituting (A-17) into the following: 

x 5 + rxg = (a o> 4 cos «j) Xg 

x 5 + rxg = awj cos o> 4 t[C sin Wgt + D sin cu^t + E cos Wgt + F cos w^t ] 



A-22a) 

(A-27) 



(A-28) 
(A 29) 



(A-30) 

(A-31) 
(A -32) 
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rx 5 + x 5 = Cawj cos c^t sin oj^t + Dac^ cos o^t sin + Eau^ 



+ Faw^ cos o^t cos cu^t 



Trigonometry Relations 

cos a^t sin o>2 1 = — [sin(a>jt + cj^t) - sinfu^t - cj^t) ] 

2 

cos Ci>jt sin o^t = — [sin(a)jt + a^t) - sin(cj^t - ct^t) ] 

2 

COS Ctfjt COS OJ^t = — [cosfct^t + (Jj^t) + cos(0)j| t - OJgt) ] 

2 

COS U)^t COS Ct^t = — [cosCu^t + UJ^t) + COS^t - o^t) ] 

2 

Wgt = (c^t + O^t) 

W 4 t = (Wjt-a^t) 

M • 1 

rx 5 + x 5 = — Cac [sin(2a> 1 t + a> 2 t) - sin(-cc> 2 t) ] 

2 

+ — Dac^ [sin(2w^t -c^ 2 t) - sin(a> 2 t) ] 

+ — Eacj^ [cos(2cj^t + a> 2 t) + cos(-a> 2 t) ] 



+ Jl Fact^ [cos(2cu^t -<x> 2 t) + cos(o> 2 t) ] 
2 



COS OJ^t COS Cl>gt 

(A -32 a) 

(A- 33a) 
(A -3 3b) 
(A-33c) 
(A-33d) 

(A- 34) 
(A-35) 



(A- 36) 



Let 

= G sin(2ct)^t + o> 2 t) + cos(2cj jt + a> 2 t) + J sin(2cj^t - u> 2 t) + K cos(2o)^t - c^gt) 
+ L sin(-a; 2 t) + M cos(-c^ 2 t) + N sin(ct> 2 t) + P cos(a> 2 t) (A- 37) 

= G(2o) 1 + c t> 2 ) cos(2a) 1 t + cj 2 0 - H 1 (2a> 1 + a> 2 ) sin^a^t + a> 2 t) 

+ J(2a> 1 -ct) 2 ) cos^a^t - a> 2 t) -K(2ct>^ -a^) sin^a^t - a> 2 t) 

- La> 2 cos(-a> 2 t) + Mcj 2 sin(^a; 2 t) + Na> 2 cos(ct> 2 t) - Pa> 2 sin(a> 2 t) (A-37a) 
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x 5 = - G(2co 1 + c t> 2 ) 2 sin(2c^2t + o^t) - + a> 2 ) 2 cos(2a>]t + u^t) 

2 9 

-J(2w 1 -cj 2 ) sin^Wjt- o; 2 t) -K(2 (j 1 -cj 2 ) cos(2cu.|t - a^t) 

2 2 2 9 

- L c ^2 sin(-aj 2 t) -Mw 2 cos(-a> 2 t) - N a) 2 sin(a> 2 t) - P u> 2 cos(a> 2 t) 

(A- 37b) 

Substitute (A-37), (A-37a) and (A-37b) into (A-36). 

2 2 
-rG(2cu 1 + cu 2 ) sin(2cu^t + Wgt) -r H 1 (2cu 1 + cu 2 ) cos^cujt +. <^ 2 t) 

2 2 

-ttJ(2cu 1 - cu 2 ) sin(2cujt - a^t) - rK(2co 1 - cu 2 ) cos(2a? 1 t - a> 2 t) 

2 2 2 2 

-tLcu 2 sin(-cu 2 t) - tM cu 2 cos(-o) 2 t) -r Nu> 2 sin(cu 2 t) - r P cu 2 cos(cu 2 t) 

+ G(2cu 1 + cu 2 ) cos(2w^t + a> 2 t) - + cu 2 ) sin(2u>]t + a> 2 t) 

+ J(2cu 1 -cu 2 ) cos (2cu 1 t - cu 2 t) - - a^) sin(2cu 1 t - a? 2 t) -Lcu 2 cos(-cu 2 t) 

+ Mcu 2 sin(-a> 2 t) + Ncu 2 cos(cu 2 t) - P a> 2 sin(cu 2 t) = 1/2 Cacu^ sin(2cu^ + w 2 t) 

- 1/2 Cacuj sin(-cu 2 t) + 1/2 Dacu^ sin(2ajjt -a^t) - 1/2 Dacu^ sin(cu 2 t) 

+ 1/2 Eao? 1 cos(2cu^t + a^t) + 1/2 Eacuj cosf-a^t) + 1/2 Facu^ cos(2cu 1 t -cu 2 t) 



+ 1/2 Fao? 1 cos(o> 2 t) 






(A- 38) 


Now equate like parts and solve for coefficients. 

9 

sin(2ojjt + a> 2 t): -tG(2oj^ + Wg) — H^(2cu^ + ^2) = 


— Caw. 
2 1 


(A- 39) 


cos (2cu^t + w<£) : 


2 

-tHj(2cu^ + ojg) +G(2 o> 1 + cu2) = 


— Eacu., 
2 1 


(A -40) 


sin (2cujt - o^t): 


-r J(2co^ - cug) 2 - K(2o)^ - cu 2 ) = 


— Daw. 
2 1 


(A-41) 


cos(2cu 1 t- o) 2 t): 


-rK(2cu 1 - a^) 2 + J(2ajj - cug) = 


— Faco., 
2 1 


(A-42) 


sin(-u> 2 t): 


-rLo)^ + Mcj 2 = - 1/2 Cacu 1 




(A -43) 
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(A -44) 



cos (- a> 2 t) : 

sin (u 2 t) : 
cos (a? 2 t) : 



o 

-rM o> 2 - Lc t> 2 
2 

-tNo) 2 -Po> 2 
2 

-rP o? 2 + No> 2 



= — Eao> 1 
2 1 

= -- Daco, 

2 1 

= — Facd* 
2 1 



H* 



o 

- 1/2 Cacu^ - rG( 2ca)^ + cu 2 ) 
(2o)^ + cu 2 ) 



(A45) 

(A-46) 



(A- 39a) 



Substitute (A- 39a) into (A-40) and solve for G: 



-r (2cu^ + cu 2 r 



- 1/2 Cac^ - t G^cu^ + a> 2 r 

(2o>^ + a> 2 ) 



+ G(2cu^+cu 2 ) =— Eacu^ 
(A-40) 



— rCao) 1 (2cu 1 + o> 2 ) + r- G(2o> 1 + + G ( 2a >i + ^ = — Eaw i 

2 2 

(2«^ + o> 2 ) + r 2 (2wj + d> 2 ) 3 = — Eawj rCaoij (2a> 1 + a> 2 ) 

I 2 2 

Eacu* — 7"Cact)| (2W| + cu 2 ) 



G = 



2[(2cj 1 + cu 2 ) + r 2 ( 2 cj 1 +cu 2 ) 3 ] 



(A-47) 



or 



G = 



G = 



acu 1 [E— rC( 2u)^ + u> 2 ) ] 
2[(2cu^ + (^2) ■*“ T (2w^ + cu 2 ) ] 

2 

l/2Eacu 1 + t H 1 (2o> 1 + c 



-r( 2u)^ + cu 2 )^ 



(2 cOj + u> 2 ) 

into (A -39) ; 

1/2 Eau)j + r H 1 (2o> 1 + u) 2 y 



(A-47) 



(A -40a) 



Substitute (A -40a) into (A- 39) and solve for H x : 

x 2 



(2cu^ + a> 2 ) 



- H 1 (2cu 1 + cu 2 ) = — Cacu^ 



(A -39) 
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1 o 1 

- — rEacOj (2co^ + o> 2 ) - r + co 2 ) - H^(2 cj^ + u^) = — Cau^ 

2 2 

p q i i 

H 1 [ (20^ + co 2 ) + r (2 cl> 1 + cj 2 ) ] = Caco 1 r Eaco 1 (2o> 1 + a> 2 ) 

2 2 



Hi = 



-Cacoj -r Eacj 1 (2cu 1 + <u 2 ) 

2 [ (2 Ct> + a> 2 ) + (2a) ^ + a> 2 ) 



(A -48) 



or 



H i = 



aa^ [- C -r E^coj + u> 2 ) ] 
2[(2co^ + co 2 ) + t^(2cj^ + a> 2 ) 3 ] 



J = 



l/2Faa>i+ T ^(2^i~^p) 



1 “ w 2' 



-t(2o>i - a ; 2 r 



(2oj 1 - a> 2 ) 

i solve for j 
1/2 Faa?i + rK(2w^ - c 



Substitute (A -42a) into (A- 41) and solve for K 

v 2 ' 



(2co x - w 2 ) 



(A -48) 



(A -42a) 



- K(2 Wl - w 2 ) = i- Dao>! (A-41) 



— — r Fact>i (2a)i - co 2 ) - t^K(2wi - a; )^ - K(2wi - w 2 ) =— Dawi 
2 2 

K[ (2co^ - cu 2 ) + r^(2cu^ - o> 2 ) 3 ] = - — Dacu^ - — r Facu-^(2co^ - cu 2 ) 

2 2 



aa^ [- D - t F( 2 cuj - cu 2 ) ] 

K. = ... , — . . ....... . " ■ ' • ■ — 

2[(2o> 1 - cu 2 ) + r 2 (2a? 1 - cj 2 ) 3 ] 



(A -49) 



K = 



- 1/2 Dao)^ - r J(2a)^ - a> 2 )^ 
(2a) 1 - cj 2 ) 



(A-41a) 
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Substitute (A-41a) into (A-42) and solve for J ; 



-t (2cu j - co 2 )^ 



- l/2Dao> 1 -t J(2cu 1 - cu 9 ) 



1 w 2' 



(2c^ 1 — cu 2 ) 



+ J(2o)j - cu 2 ) = — Facu. 



— r Dacu^(2cu 1 - o? 2 ) + t 2 J (2cUj - cu 2 )^ + J(2cu 1 - o> 2 ) = — Faco^ 



2 3 1 1 

J[(2cu^ - cu 2 ) + t (2co^-o> 2 ) ] = — Fac rDacUj (2cu^ - co 2 ) 

2 2 



J = 



acu^ [F - tD(2cUj - a> 2 ) ] 

2[ (2cu^ - o> 2 ) + t 2 (2cu 1 - o> 2 )^ ] 



M 



l/^CacUj + tL c^2 



Substitute (A- 43a) into (A- 44) and solve for L ; 



-T CUo 



- l/2Cacu 1 + rL ( 



- Lcu 0 = — Eacu 1 
z 2 1 



,2 T 3 



1 



— r Cacu 1 cu 2 -r La> 2 - Lo> 2 = — Eaco 



2 3 1 1 

L(ct> 2 + 7 cu 2 ) = - — Eacu^ + — r CacUjO^ 

2 2 



L = 



awjJ- E + r Cco 2 ] 



2[ Ct>2 ^ k>2 1 



L = 



- l^EacUj - cu 2 



0) o 



(A-42) 



(A- 50) 



(A-43a) 



(A -44) 



(A- 51) 



(A- 44a) 
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Substitute (A-44a) into (A-43) and solve for M 



-T CU 0 



- 1/2 Eacu^ - r M o>2 



Mcu 0 = CauK 



1 2 3 1 

— TEacu^ct>2 + r Mo> 2 + Ma>2 Cacu^ 

2 3 1 1 

M(o>2 + r cu^) = Cacuj rEacu^^ 



M = 



ac (- C - r Eo) 2 ) 
2(ct>2 + ^2^ 



l^DaWj- rN c^2 



Substitute (A-45a) into (A-46) and solve for N: 



1/2 Dac^-rN 



+ Ncu 0 = — Facu., 
1 2 1 



1 2 3 1 

rDaco^cu2+ t Ncug + No> 2 = — Faco^ 



N(a>2 + r 2 cu|) = — FacUj + — r Dacu^ 
2 2 



N = 



au> 1 (F+ tDcu 2 ) 
2 3 

2(o>2 + t a> 2 ) 



N 



1/2 Faa^ + r P 



(A-43) 



(A- 52) 



(A- 4 5a) 



(A-46) 



(A- 53) 



(A-46a) 
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Substitute (A -46a) into (A -4 5) and solve for P : 



1/2 Fac Uj + t P cug 



-Pcu 2 = 



1 



- — Dacu 



1 



1 2 3 1 

r Faa> 1 cu 2 - r Pw 2 - Pa> 2 = Dacjj 



2 3 1 1 

P(cu 2 + r cu 2 ) = — Dacu 1 r Fao^cug 



acu^ (D - r Fcu 2 ) 



2(ct> 2 + T 2 Ct> 2 ) 



X = Xj + Xg + Xg + 



^m-l + *2 



m 



x^ = A sin cu 2 t + B cos u> 2 t 



Xg = C sin cu^t + D sin cu^t + E cos cu^t + F cos cu^t 



(A-45) 



(A- 54) 



(A- 3) 
(A-7) 

(A- 17) 



x & = G sin (2^^ + a> 2 t) + Hj cos (20^ + cu 2 t) + J sin (2a?^t - cu 2 t) + K cos - cu 2 1 
+ L sin(-ct> 2 t) + Mcos(-a) 2 t) + Nsin(co 2 t) + Pcos(o> 2 t) (A- 37) 



Terms beyond x^ may be neglected if the order of magnitude of a is much 

3 

smaller than unity (maximum a on the MIT angular shaker is 2 x 10" 

radian). If x^ is defined to be of the order of magnitude of one, x^ is of the 

2 

order of magnitude of a, and x^ is of the order of magnitude of a . Hence 
each successive term in the series of x is smaller than its predecessor in 
accordance with the value of a. 

For further discussion on convergence and a process for iteration, 
ref. 10 is recommended. 



73 



If only the terms associated with Xp x^, and Xg are considered, it is 
readily apparent that the time average of drift for the sinusoids is zero. Close 
scrutiny, however, shows drift possibilities as and a> 2 approach one another 
in value and also when harmonics are present. 

The specific case where = u> 2 has been investigated in ref. 10. The 
general case which is the subject of this chapter degenerates to the specific 
case when = o> 2 . Under such circumstances D sin t and F cos a^t 
require L' Hospital's rule for proper evaluation, but eventually yield a drift 
term which is a direct function of time. 

The other case of interest concerns harmonics. The series terms asso- 
ciated with Xg require L' Hospital's rule for evaluation when 2u>^ = a> 2 . Under 
such a condition drift terms appear which are direct functions of time. 

An example of drift caused by the existence of harmonics such as when 
2o>^ = a> 2 is shown for the term J sin(2o>jt - a> 2 t) . 

If 20^-* u> 2 , then - a> 2 ) -> 0 (A-55) 

au)., [F - tD (2cu 1 - cu 9 ) ] 

J sin (2o>^t - u> 2 t) = = = sin (2o> 1 t - a> 2 t) (A- 56) 

2 (2w 1 - w 2 ) [1 + r 2 {2w x - w 2 ) 2 ] 

Let (2a>jt-u> 2 t) = (2a)j-a> 2 )t = yt (A- 57) 

lim J sin yt = - 9 - apply L' Hospital's Rule (A- 58) 

y+0 0 



— [acu^F- rDy) sin yt] 
dy 

— [2y(l + r 2 y 2 ) ] 

dy 



ac 

~ 2 ~ 



- t D sin yt + (F - rDy)(t cos yt) 
(1 + t 2 y 2 ) + y (2 r 2 y) 



(A -59) 



74 



lim J sin yt 

y-»o 



aw l (F)(t) 



Faco^t 



(A -60) 



2 1 



2 



If terms for x^, Xg, etc. were evaluated, drift terms would appear for higher 
harmonics. As the series has an infinite number of terms, it seems apparent 
that any two given input frequencies might create some type of harmonic drift , 
but the magnitude of "a" associated with such drift may make the consideration 
of the value a folly when compared to the manufacturing tolerances used when 
constructing gyros. 

The case of negligible elastic restraint shows the possibility of drift terms. 
However, for such a gyroscope, the shortened series representing drift is valid 
only when values of a, b, and c are small compared to unity. Furthermore, 
man can only produce a gyro which is not perfect. Hence terms in the series 
which are of less consequence than the manufacturing tolerances will have only 
a theoretical importance and should not enter practical considerations. 

A. 2 Series Solution for Finite Elastic Restraint Between 
the Float and Case About the Output Axis (k 4 0) 

The general equation is : 

(k - aa> i cos co^t) x + x + r x = bcog cos c^t - ca^ sin + kb sin o^t 



In accordance with the series established in Section 5. 6, the interest lies with 
the following terms : 



+ kc cos co <£. (A-l) 




(A-3) 



where: 



kx^ + Xj + TXj = ba »2 cos cogt - eu >2 sin a^t + kb sin cogt + kc cos cogt 



(A-61) 
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(A-62) 



• • • 

kx^ + x 3 + rx 3 = (ac*^ cos o>^t) x x 

kX[5 + x 5 + rx 5 = ' aa, i cos ^1^ x 3 

• •* 

(k-aco 1 cos cjjtJx + x + rx = bcu 2 cos co 2 t - ca> 2 sin cu 2 t + kb sin cu 2 t 

+ kc cos w 2 t 

Let: 

x = Xj + X 2 

where: 

• •• 

k x i + Xj + tXj = ba; 2 cosa) 2 t - ca> 2 sincu 2 t + kbsina> 2 t + kccosa> 2 t 
Substitute (A -64) and (A-65) into (A-l). 

t • •• M 

(k- aco^ cos cu^t) (x^ + x^ + (x^ + x 2 ) + r (x^ + x^ = ba> 2 cos a> 2 t - cci> 2 

+ kbsincu 2 t + kccoscu 2 t 
or 

• •• 

(k - aa^ cos cu^t) X 2 + x 2 + tx 2 = (acu^ cos cu^t) x ^ 

Try a solution for (A-64) of : 

x^ = A sincu 2 t + B cos cu 2 t 

x^ = A cu 2 cos cu 2 t - B ct> 2 sincu 2 t 

*•2 2 

x^ = Ao> 2 sincu 2 t - Bcu 2 coscu 2 t 

Substitute (A-67a) and (A-67b) into (A-65). 

2 

k A sin a> 2 t + k B cos cu 2 t + A a> 2 cos a> 2 t - B a> 2 sin a> 2 t - t A a> 2 sin a> 2 t 
o 

- TBa> 2 cosa> 2 t = bwgcoswgt - ca> 2 sinw 2 t + kbsinci> 2 t + kccosa> 2 t 
Equate coefficients of like terms. 

2 

sin • k A - - r ^ w 2 = - CU> 2 + ^ 

2 

cos cl >gt : kB + Aco 2 - tBo> 2 = bt> 2 + kc 



(A-63) 

(A-l) 

(A-64) 

(A-65) 

, sin Wgt 
(A -66) 

(A-66a) 

(A-67) 

(A-67a) 

(A-67b) 

(A-65) 

(A -68) 
(A -69) 
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(A-68a) 



2 

Afk-rcog) = Bo> 2 -co >2 + kb 



tB o)« -kB + bo> 0 + kc . ^ , 

A = ? 2 . tBu , 2 _ + b + kc_ 



Wo 



v 2 ^2 
Substitute (A- 69a) into (A-68a) and multiply through by ct> 2 • 

kTBcu^ -k 2 B + kba>2 + k 2 C - T 2 Ba>2 + rkBc^ - Tbco 2 ” ^kccu^ 

2 2 

= Bco 2 - cu>2 + kbcu 2 

B[ Tcjg - k 2 - t 2 o>2 + T.kcUg - ^ ] = T bct >2 - k 2 c + r kc cu 2 - co )2 



Ct)n 



B 




1 + 




- r a>2 c I — + 



_b 







B 



Ak — At a>2 + ccug - kb 



Ak _ Arct> 2 + c — 

w 2 w 2 



Substitute (A-68b) into (A-69a) and solve for A : 



A = 



( r a>2 - k) 



Ak » „ kb 

- At oo 2 + c 



+ b + 



kc 

Wo 



(A -69a) 



(A- 70) 



k 

w 2 



(A-71) 



(A- 68 b) 



(A- 69) 
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2 2 

°2 T ' "2 ' 



Ak 



kc kb 



W 9 ,2 
2 w 2 



A 1 -r k + r 2 o >2 + | 2 -r k 


k 2 b 

- Tea, -Tkb + + 

z 2 
w 2 




b I 1 + f I + T UJnC 

A = |_ K/ J 


c \ w 2 lj 




1 + 


to ^ 


| - T w 2 


2 




(k-acu^ cos cu^tjxg + x 2 + tx 2 = (ao>i cos cu^t) x^ 



Let: 

where: 



x 2 = x 3 + x 4 



kxg + x 3 +rx 3 = (ao^ cos co^t) Xj 
Substitute (A -7 3) and (A- 74) into (A -6 6a). 

(k - aw 1 cos w 1 t)(x g + x 4 ) + (x g + x 4 ) + r (Xg + x 4 ) = (awj cos a^t) x 4 (A-66a) 

or 



kc 

CUn 



(A- 72) 

(A -6 6a) 
(A-73) 

(A- 74) 



(k - aa>i cos o^t) + x^ + r = (aco^ cos co^t) Xg (A-66a) 

A solution is desired for the following equation: 
kXg + Xg+TXg = acu 1 A CO s O) 1 1 sin cu 2 1 + a co 1 B cos co 1 1 cos co 2 t (A-74a) 

Let: 

x^ = D sinojjt sino^t + E sincu^t cos cugt 

+ Fcos cu^t sin cugt + Geos cujt cos cugt (A-75) 
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Xg = cos ojjt sina^t + sin a^t cos a^t 
+ Eo>^ cos cu^t cos o^t - Ko >2 sin cu^t sin o^t 

- F^ 1 sin cjjt sin o^t + Fu> 2 cos cUjt cos c^t 

- sin cj^tcoscjgt - Gcjg cos Wjt sincjgt (A-75a) 

•• 2 

Xg = -Dcu^ sin co^t sin cUgt + Da^a^ cos cos a^t 

2 

+Du)^2 cos cu^tcos a^t-DcOg sincu^t sina^t 

2 

-Ea^ sina^tcoscogt - Ea^cug cos co^t sina^t 

2 

-EcUiO^cosa^tsina^t - Ea^ sin a^t cos a^t 
2 

-Fajj cos cu^t sina^t - Fa^a^ sin o^tcosc^t 

2 

-Fcu^a^ sin cu^t cos c^t " FcUg cos cu^t sin a^t 

2 

-Gc^ cos cu-^t cos a^t + Ga^cUg sincu^t sin a^t 

2 

+ Ga^c^ sincu^t sina^t - G^ cos cu^t cos 0^2^ (A- 7 5b) 

Now substitute (A-75), (A-75a) and (A-75b) into (A-74a) and solve for D, E, 

F, and G. 

(kD sin o^t sin Cc>2t + kE sin cu^t cos a^t + kF cos cu^t sin a^t + kG cos cu^t cos a^t) 

+ (Dcu^cos cu^t sin u^t + Da>2 sin cos a^t + Ecu^ cos co^t cos a>2t - Ea>2 sin cu^t sin a^t 
-Fa>^ sin cu^t sin cugt + FcUg cos cos a^t - Gcu ^ sin cu^t cos a^t -Ga^ cos a^t sin a^t) 
2 

+ (-rDcu^ sincu^tsin^t + tDcu^o >2 coscu^tcosc^gt + rDajjWg cos cu^t cos cugt 

2 2 

-rDcu 2 sincu^t sin cu 2 t - tE sincu^t cos a^t - r E cos cu^t sin u^t 

2 2 

-TEa>^cu 2 coscu^tsincugt - rE cu 2 sin cu^t cos cu 2 t - r F cos a^t sin a^t 

2 

-rFcu^cUg sina>^tcos a^t - rFa^a^ sina>jtcos a^t - tF<a >2 cos cu^t sina^t 

2 

— rG cos cu^t cos a^t + rGcojOJg sincu^t sina^t + tGcu^o^ sin cujt sincugt 

2 

-rG cug cos cujt cos a> 2 t) = Aacu^ cos o^tsino^t + Bac*^ cos a^t cos a^t (A-76) 
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Now equate like terms of the various trigonometric function combinations 
and solve for D, E, F, and G. (Note: A and B are known) 

(sin co^t sin o^t) 

2 2 

kD - Eo>2 " Fa, l " ^1 ” + t ^ Cx) i U) 2 + r ^ co x co 2 = ^ 

kD - Eco 2 - Fo)j - D(t + r o> 2 ) + 2r Gco^o^ = 0 
D(k- r - r cu 2 ) - Eo? 2 - Fcu^ + G(2 rcu^o^) =*.,0 (A-77) 

(sincu^t cos co 2 t) 

2 2 

kE + Dcu 2 - Gcu^ - rE - tE a> 2 - rFajjWg - rFco^cu 2 = 0 

2 p 

Dco 2 + E(k- r - r cu 2 ) - F(2 tcu^co 2 ) - Gcu^ = 0 (A-78) 

(cos cu^t sinco 2 t) 

2 2 

kF + Dco^ - Gcu 2 - t Eco^cu 2 - rEco^cu 2 - rF - rFcu 2 = Aaco^ 

2 2 

Dco^ - E(2 t co^ct> 2 ) + F(k - r - r cu 2 ) - Go> 2 = AacUj (A-79) 

(cos cu^t cos o> 2 t) 

2 2 

kG + EcUj + Fcu 2 + tDcOjCi> 2 + tDcu^cu 2 - rG - rG o> 2 = BacUj 



2 2 

D(2to)^o; 2 ) + Ecu^ + Fco 2 + G(k - r - r cu 2 ) = Bacu^ (A-80) 

9 2 

D(k - Td)J-T o> 2 ) - Ecu 2 - Fcu 1 + G(2 r cu^) = 0 (A-77) 

Dcu 2 + E(k — r — r Ci> 2 ) - F(2 r cu^cOg) ” = 0 (A-78) 

Dcuj + E(2 rcu^cu 2 ) + F(k - r - r cu 2 ) - Gcu 2 = Aacu^ (A-79) 

D(2 r cu 1 co 2 ) + Eco 1 + Fo> 2 + G(k - r ) = Bao^ (A-80) 

Solve for D, E, F, and G by determinants. 



(A-81) 



D = 



0 


-"2 


- W 1 


27a;^a;2 


0 


(k- ro)j - rwg) 


-270)^0^2 


- w i 


Aacoj 


27 0)^0^2 


(k- TWj - TWg) 


-a» 2 


Bao^ 


w i 


w 2 


(k - 7 - 7 a>2 ] 


to) 2 - ro)2) 


~ W 2 


- W 1 


2 7 a)^c^2 


w 2 


(k-T^-TWg) 


- 270)^2 


- w i 


W 1 


27(0^2 


(k- TO)^- TWg) 


- w 2 


2 TOJ^U)^ 


w i 


w 2 


(k- Ta)j - TWg 








(A-82) 


denominator of 


:he determinant be known as A or (A-82) . 




(k - TCJj - TWg) 


- 2tu>iu)2 


- W 1 


TWj- TUj) 


2tw^W2 


(k- To )^- TWg) 


~ w 2 




W 1 


W 2 


(k-rw^-TWj) 




w 2 


- 2705^(^2 


- W 1 


+ o> 2 


W 1 


(k-TO)j - TOJg) 


-u)2 




2 TojjWg 


w 2 


(k- Tia\- T0}\) 




w 2 


(k - t Wj - t a) 2 ) 


~ W 1 


- W 1 


W 1 


2 7(0^2 


-a> 2 




2 t oj^cl>2 


W 1 


(k- rwj - raig) 




0) 2 


(k - T CO ^ - TWg) 


- 27 OJjWg 


-2 to)^2 


W 1 


2 rajjWg 


(k- 7 0^ - 7 0^) 




2 TW^Wg 


a»i 


w 2 



(A-82) 



81 



A = (k - TMj - TWg) ^(k - TOJj - TWg) + (k- TO)j - TWg) - 2 TWj U>2 

+ w^k-TWj-ruj) + 2 tWjU> 2 + 47-2 w i w 2 (k " TW 1 " rw 2 ) ] 

[ „ 2 2 x 2 3 o 2 „ 2 2-. 

u> 2 (k- tWj - tu> 2 ) + w 2 + 2 r oij w 2 (k - r - r« 2 ) 

2232 2 2 2 1 

+ 4 r 0*1 a>2 - o^2 + 2 r co^ cog (k - - ro^ ) I 

F 2 22 2 2 2 ^ 

+ (-k>j) 2 ro> 1 co 2 (k- ro)j -tcj 2 ) + ^1^2 - ( k - T - r ^ 2 ^ 

rt 2/, 2 2* 3 .2321 

- 2 co 2 (k — r co ^ — t c^ 2 ) - Wj + 4 r a> 1 co 2 

[ 3 2 2 2 2 

2 ra^cog - a^a^k - TWj- tu> 2 ) ~ a^a^k- r - r co 2 ) 

2 ~i 

+ 2 r (k - r - r cj 2 ) - 2 r w 2 + 8 r ^ cu^ co 2 1 (A-82) 

4 2 
A = (k — t coj — r Ct> 2 ) + (2ct)^ + 2 ^ 2 ) (k - r - r a?^) 

+ 8 r wg (k - r coj - r u 2 ) - 2 a> 2 a> 2 + a)* + o)\ - 16 f 2 a>* co* (A-82a) 



Determinant Form of Coefficients D, E, F, and G. 




(A-83) 



82 



0 



2 to>iU>2 



F = 



(k - T^ - T COg) 



-C 0, 



Wi 



2 tco^cj2 



(k - T GJ 2 — TO) 2 ) 



2 tcj^cj2 



Aaa>i 

Bac^t 



-Wi 



-C0 o 



2 2 x 

(k- TC 0 1 - TWg) 



(A-84) 



G = 




AD = Aawj j^rWjW^k-rw^-T wg) + 4 t 2 (j|oi 2 2 + wJ-WjW 2 
9 9 2 2 2 \ 

+ 2 T 0^0*2 (k- TWj -r o> 2 ) + o? 1 (k - r - r a> 2 ) J 

2 2 2 

-Baw 1 j^-2r WjWg + 2 tw 2 w 2 + 2t Wj&^k-r Wj - tu> 2 ) 

+ 8twj(d 2 - w 1 w 2 (k- m] - t w 2 ) -w 1 w 2 (k- r u> 2 - t w|)] 

(A-81) 

AD = Aawj £ 4 t C 0j(t> 2 (k — r « 2 -t w 2 ) + 4 t 2 o>x« 2 + w x - w i w 2 
+ Wj (k - r w 2 - r « 2 ) j + Baw^ ^2 tWjW 2 - 2 r w 2 



- 2 r WjW 2 (k - t w 2 - t « 2 ) - 8 r to 2 w 2 + 2w^w 2 (k - r 



2 2 J 

'1-WJ 



(A-81) 



83 



[ 2 2 3 3 

- 2 re* ^2 (k - r o> 1 - r ct> 2 ) + 2 t cjg + 2 rco^g 

333 22 2 2 1 

+ 8r cog + o^o) 2 (k- rco| - tco 2 ) + o^a^k- rco^ - rcop I 

[ 2 2 2 3 2 2 2 

2 (k - tcj^ - rcjg ) + + 2 Tco^cOg (k - rco^ - rct> 2 ) 

+ 4t 2 w 2 u| + w^k-rwj-T^) + (A-83) 

2 

AE = Aac^ £ 2 rcu^c^k- T c^ - rci> 2 ) - 2 rc^ c^ 2 - 2 - 8 co^ co 2 

2 2 1 r 2 223 

- 2o>^ct?2(k - rWj - rcj 2 y + Bacu^ 4 tco^o^ (k - rcu^ - r ci> 2 ) + coj 

2 

+ 4 t 2 cog + ^(k - to)^ - rcog) + c^a) 2 J (A-83) 

[ 2 2 ^ 2 2 2 2 

(k - tci>^ - rct> 2 ) + 2 cog + 2 ro) 1 cog 

.222 /, 2 2 x 2/. 2 2x 2,. 2 2 X 1 

- 4 T Ct)j CUg (k - TCOj - T(x) 2 ) + COj (k - TWj - TCJg) + Ct> 2 (k - TCO^ - TCOg) 

9 J 

[ 2 2 ^ 2 2 2 3 

— Ct> 2 (k — t — t Ct> 2 ) + co 2 + 4 r o> 2 

- 2 rct)^ cjgCk- - rcjg) + 2 ro)^ a> 2 (k- ra^ - tcj^) - co 2 j 



[ o 0 ,5 , 22 / Zyizzz Z\ 

(k-rw 2 -rw 2 ) + 4 twj a> 2 + (c^ - 4 r WjO^ + Wg) 



.2 2,2 



(A-84) 
. 2 , 



(k- ra ) 2 - t C t> 2 ) 

.2 2 3 3 

■ 4 r co 2 + a> 2 



[ 2 2 2 o 

co 2 (k - rct)^ - ra> 2 ) - co 2 



(A-84) 



[ 2 2 ^ 2 2 3 2 

a> 2 (k - to*! - ro) 2 ) +4r a> 2 - 

+ 2 rcoj a? 2 (k - - tcu 2 ) + 2 ra^ a> 2 (k - ra^ - rco^ ) + cu 2 1 

3 J 

+ Bact^ £ (k - ro)^ - ro> 2 ) +2 tc^ o> 2 - 2 ro)^ k? 2 + (k - t < jo ^ - ro) 2 ) 

+ 4 w| o> 2 (k - ro?^ - to? 2 ) + o> 2 (k - rco^ - ra> 2 )J ( A-85) 



84 



AG = Aao> 



i[ w i w : 



2 - w 2 (k- 70)2-70)2) 



.2 2 3 
4 r c^2 



2 2231 T 2 2 ^ 

- 4 rcjj a) 2 (k- rco 1 - r^) - c^2 + Baa, i (k - - rcug) 

+ (cu^ + 4 r 2 o>2 + u\) (k - r - rcu^j (A-85) 



(k - aa)j cos cu^t) x 4 + *4 + rx^ = (aco^ cos co^t) x^ (A-66a) 

x 4 = x 5 + x 0 (A-86) 

x 3 = D sin a>jt sin + E sin cu^t cos (o^l 

+ F cos ajjt sin u>2 1 + G cbs a>^t cos o^t (A - 7 5) 



where. + + rX g = (acu^ cos cu^t) x^ (A-63) 

Substitute (A-86) and (A-75) into (A-66a). 

• • •• •• 

(k - ao>i cos cu^t) (Xg + Xg) + (x^ + Xg) + r(Xg + Xg) 

= Dacuj sin cu^t cos u^t sin u> 2 t + Eao> 1 sin co^t cos co^t cos a> 2 t 
2 2 

+ Facu^ cos co^t sin c o<£ + Gacu^ cos cujt cos cu 2 t (A-87) 



Try a solution for x^ where: 



kx 5 +x 5 + tx 5 = 



Daco.. Eacu.. 

i sin 2c«)jt sin cu 2 t + ± sin 2o? 1 t cos o^t 

2 2 



Fao)j 

2 

Gao> 



Facu- 



sin 0) o t + 



cos 2(o A sin cu 0 t 



Gacu 1 



— cos cu«t + cos 2co A cos co 0 t 

2 2 2 1 2 . 



(A -88) 



Let 

x 5 = s * n ^1* sin w 2^ + 1 sin cos 
+ J cos 2co^t sin o^t + L cos 2 co^t cos cogt 

+ M sin o> 2 t + N cos C0 2 t (A-89) 



85 



Xg could be written in an equivalent form with certain terms containing 
(2o)^t - a?2 1) by use of trigonometry relations. Thus a time dependent 
drift term is seen to appear when the harmonic situation = cjg 
exists. 

Xg = 2Hj cos2a>jt sin co^t + cu 2 sin 2o>jt cos a^t 
+ 21 cos2cu^tcoscu 2 t * I cug sin2cu^t sinc^t 

- 2 J u> i sin 2cu ^t sin + J cug cos2cu^t cos a^t 

- 2L sin2cu^tcos cjgt - L co 2 cos2cu^t sincugt 

+ Mcugcoscugt - NcUgSincjgt (A-89a) 

Combining terms of x^, one obtains: 

x & = (2H 1 cu 1 - La^) cos 2co l t sin C0 2 t + ^l co 2 “ 2Lcj i) sin 2cu 1 t cos c^ 2 t 
+ (2Icu^ + JcUg) cos 2cu^tcos cugt - (Ia^ + 2Jcu^) sin2cu^tsin a^t 
+ (Mcu 2 ) cos cjgt - (Ncu 2 ) sin cu 2 t (A-89a) 

x^ = (2HjW^ - La> 2 ) (“2 cu^) sin 2cu^t sin a> 2 t + (2H^cu^ - Lcjg) ^ 2 ) cos 2cj i^ cos 
+ (H^cj 2 - 2Lcu^) (2cup cos 2cu^t cos cu 2 t + (H^cj 2 - 2Lcj^)(-cu 2 ) sin 2cu^t sin cu 2 t 
+ (2Icu^ + Jcu 2 )(-2cj^) sin2cu^t cos cu 2 t + (2Icu^ + Jcu 2 )(-cu 2 ) cos 2cu^t sincu 2 t 

- (Ia>2 + 2Jcj^)(2cu^) cos2cu^t sincu 2 t - (Icu 2 + 2Jcu^)(cu 2 ) sin2cu^tcos cu 2 t 

- (Mcu 2 ) sincu 2 t - (N cu 2 ) cos cu 2 t (A-89b) 

Combining terms of x^, one obtains: 

x 5 = (2 Lcu^^ 2 “ 4Hj cu 2 - H 1 w 2 + 2La> 1 cu 2 ) sin 2cu 1 t sin cu 2 t 
2 2 

+ (2H^cu^cu 2 — L cu 2 + 2H 1 cj 1 cj 2 - 4Lcu 1 ) cos 2 ( 0 ^ cos a^t 

2 2 

+ ( — 41 y ^ — 2 JcUjCU 2 — 1 ^2 — 2J(02^2^ sin 2cu^t cos CcJ 2 b 
+ (- 210)^2 - Jcjg - 210)^2 ~ 4Jcu^) cos 2cu 1 t sin cugt 
+ ( - M cu 2 ) sin cugt + ( - N cu 2 ) cos cUgt (A-89b) 



86 



Now substitute (A-89), (A-89a) ? and (A-89b) into (A-88) and solve for 
coefficients H^I, J, L, M, and N. 

j^kH^ sin2o)^t sin o^t + kl sin 2a) ^t cos a^t + kJcos2cUjt sinu> 2 t 
+ kL cos2o? 1 t cos cjgt + kMsincUgt + kNcosa> 2 tJ + 

£(2 H^i - Lc^ 2 ) cos 2a^t sin cu 2 t + (H^a) 2 - 2La>^) sin2co^t cos cjgt 
+ (21 cuj + Jcu 2 ) cos ^ ^t cos cjgt - (I cjg + 2J cu^) sin 2a^^t sin a> 2 t 

+ (M cu 2 ) cos aj 2 t - (Na> 2 ) sincu 2 t J + 

[ 2 2 

r (2L a > ^2 - 4Hj cu 2 + 2L cj^c sin 2^ 1 t sin cu 2 t 

2 2 

+ t(2Hj a^a> 2 - Lct? 2 + 2Hj cj^cj 2 - 4L ) cos 2a>^t cos a> 2 t 
2 2 

+ r( - 41 a>£ - 2J ct>i<^ 2 - I a> 2 - 2J a^cu 2 ) sin 2a)^tcos co 2 t 
2 2 

+ r( - 2 1 cu^a> 2 - J c^ 2 - 2 Icj^cj 2 - 4J cu^) cos 2a>^t sina> 2 t 

2 2 1 
+ r (- M a> 2 ) sin a? 2 t + r ( - Na> 2 ) cos cu 2 t J = 




Now equate like terms of the various trigonometric functions and solve 
for Hp I, J, L, M, and N. (Note: A, B, C, D y E ? and F are known) 



sin 2cjjt sin <u 2 t 

2 2 Dac o •* 

kHj - (Icu 2 + 2J cu^) + t(4Lcu^cu 2 - 4H^ o> 2 = 

» n Daa) t 

(k - 4 - rcu 2 ) + (-cu 2 ) I + (- 2cup J + (4 rcu^cug) L = 



(A-91) 
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